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INTRODUCCION 4 


MN 


ANTECEDENTES HISTORICOS 


L interés por el procesado digital de imágenes se divide 
en dos áreas principales de aplicación: mejora de imáge- 
nes (para su mejor interpretación por el hombre) y pro- 
cesado de imágenes (para su interpretación automática 
por medio de una máquina). Las primeras técnicas de tra- 
tamiento de imágenes se remontan al principio de los años 
veinte, cuando se trataba de mejorar las fotografías de pe- 
riódicos enviadas por cable submarino entre Londres y 
2 Nueva York. El uso de imágenes en forma digital se re- 
monta a 1921, con la introducción del primer sistema de transmisión de 
fotografías por cable submarino desde Londres a Nueva York. La introduc- 
ción del sistema Barlane de transmisión fotográfica por cable redujo el 
tiempo empleado en transportar una fotografía a través del Atlántico de 
una semana a menos de tres horas. Los primeros sistemas Barlane eran ca- 
paces de codificar imágenes en cinco niveles distintos de intensidad y esta 
capacidad fue aumentada en 1929 a quince niveles. Durante este período 
también se mejoró notablemente el proceso de reproducción, que pasó de 
ser la impresora telegráfica a un sistema con placa fotográfica impresio- 
nada con rayos de luz modulados por la imagen codificada. 


Aunque continuaron las mejoras en los métodos de procesado durante 
los siguientes treinta y cinco años, la gran potencia de cálculo de los com- 
putadores, derivada de la progresiva evolución de las tecnologías VLSI y 
el programa espacial fueron los que dieron gran impulso a los conceptos 
de procesado de imagen. El tratamiento digital de imágenes entendido 
como la aplicación de procedimientos matemáticos tendentes a corregir 
las fotografías y mejorar su interpretabilidad no se inició hasta 1963, cuan- 
do en los laboratorios de investigación de IBM se empezó a procesar me- 
diante ordenador las imágenes obtenidas por el satélite Tiros. A partir de 
ese momento, el tratamiento digital de imágenes se ha convertido en un 
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área de las ciencias aplicadas que viene experimentando un crecimiento 
acelerado durante los últimos veinte años. El uso de técnicas de mejora 
de imagen por ordenador comenzó en el Jet Propulsion Laboratory (Pasa- 
dena, California) en 1964, cuando fotografías de la Luna transmitidas por 
el Ranger 7 eran procesadas por un computador para corregir varios tipos 
de distorsión de la imagen inherente a la cámara de TV de a bordo. Estas 
técnicas sirvieron como base para posteriores métodos usados en la me- 
jora y restauración de imágenes para programas tan familiares como las 
misiones Surveyor a la Luna, la serie de misiones Mariner a Marte y los vue- 
los tripulados Apollo a la Luna. Este desarrollo, inicialmente, se efectuó 
con objeto de tratar las imágenes suministradas por las diversas misiones 
espaciales, en especial las de los satélites dedicados a la teledetección de 
recursos naturales e impulsada por la utilización de computadores cada vez 
más potentes, tuvo menor incidencia en otras áreas, si bien en algunos ca- 
sos, como el de la microscopía electrónica, estas técnicas se han emplea- 
do desde mediados de la década de los sesenta. 


Desde 1964 hasta ahora, el campo del procesado de imágenes ha expe- 
rimentado un vigoroso crecimiento. Además del programa espacial, las 
técnicas más usadas hoy se dedican a una gran variedad de campos como: 
medicina (imágenes biomédicas y rayos X), geografía y meteorología (aná- 
lisis de imágenes aéreas y de satélite), arqueología, astronomía, biología, 
medicina nuclear, defensa o industria. Todas ellas no sólo en cuanto a la 
mejora de imágenes para interpretación humana, sino también en lo refe- 
rente a su interpretación por una máquina en campos como: reconoci- 
miento de taracteres, robots industriales para montaje e inspección, reco- 
nocimiento de objetivos militares, procesado de huellas dactilares, análi- 
sis de muestras de sangre y radiográficas, predicción del tiempo atmosfé- 
rico y de las cosechas, acústica, comunicaciones, geofísica, procesamien- 
to de voz, microscopía, teledetección, etc. 

Por otra parte, se ha producido un gran acercamiento entre las técni- 
cas de procesado digital de imágenes y la de gráficos por ordenador (com- 
puter graphics), hasta el punto de que ambas se solapan en el marco de cier- 
tas aplicaciones. 


El tratamiento digital de imágenes abarca una gran variedad de técni- 
cas y herramientas matemáticas que se han desarrollado, si bien cabe se- 
ñalar que en los últimos años se han comenzado a desarrollar bases cien- 
tíficas en lo que antes era un simple conjunto de técnicas. 


El tratamiento digital de imágenes abarca una gran variedad de técni- 
cas y herramientas matemáticas que pueden ser clasificadas en los siguien- 
tes puntos: 


1. Digitalización y cuantificación, que consiste en efectuar el paso de 
imágenes continuas a discretas. 
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2. Operaciones de preprocesado, consistentes en la realización de de- 
terminadas operaciones tendentes a normalizar los datos de entrada. 


3. Mejora o realce de imágenes, consistente en efectuar algunas ope- 
raciones tendentes a mejorar la apariencia de los datos de entrada. 


4. Restauración de imágenes, consistente en la recuperación de la 
imagen original, en el supuesto de que los datos de entrada se encuentran 
afectados de algún tipo de degradación, como, por ejemplo, desenfoque, 
turbulencia atmosférica, movimiento, etc. 


5. Codificación de imágenes, consistente en la compresión de la in- 
formación con objeto de reducir los costes de transmisión de las imáge- 
nes, así como conseguir un almacenamiento más eficiente, como sucede, 
por ejemplo, en el caso de las imágenes médicas, en ingeniería, etc. 


6. Análisis de imágenes, consistente en la obtención de descripciones 
numéricas o simbólicas de la imagen de partida. Las descripciones de una 
imagen se refieren generalmente a las partes constituyentes de dicha ima- 
gen. Es por ello que, el primer paso en todo proceso de análisis de imáge- 
nes lo constituye la denominada segmentación en regiones con significa- 
do propio. 


7. Reconstrucción de imágenes, consistente en la obtención de infor- 
mación acerca de un objeto a partir de sus proyecciones, como, por ejem- 
plo, la realizada en tomografía, en donde se obtienen secciones transver- 
sales del cuerpo humano, o los métodos basados en la resonancia magné- 
tica nuclear. 


Estas técnicas adoptan en la práctica formas y modos de operación muy 
diversos en función del tipo de imágenes a que se aplican y del tipo de in- 
formación que de ellas se desea extraer. 


REPRESENTACIÓN 
DIGITAL DE IMAGENES 


Antes de comenzar con la descripción de algunas de las técnicas del 
procesado digital de imágenes, es preciso definir qué se entiende por una 
imagen digital, a fin de poder apreciar las posibilidades de estas técnicas, 
así como comprender la problemática que entraña su utilización. 


En el ámbito del tratamiento de imágenes, el término imagen se refiere 
a una función bidimensional de la intensidad luminosa f(x,y), donde x e y 
denotan coordenadas espaciales y el valor de f en cada punto (x,y) es pro- 
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porcional al brillo, intensidad o nivel de gris de la imagen en ese punto. 
A veces es útil visualizar una función imagen en perspectiva tridimensio- 
nal, siendo la intensidad el tercer eje. 


Una imagen digital f (x,y) se obtiene a partir de la radiancia de un ob- 
jeto que ha sido discretizada en ambas coordenadas espaciales y en inten- 
sidad. Podemos considerar una imagen digital como una matriz cuyos ín- 
dices de fila y columna identifican un punto de la imagen y el correspon- 
diente elemento de la matriz representa el nivel de gris de ese punto. Los 
elementos de dicha matriz digital se llaman elementos de imagen o pixels 
(o pels), como abreviatura de «picture elements». 


Aunque el tamaño de una imagen digital varía con la aplicación, hay 
grandes ventajas en elegir matrices cuadradas cuyos tamaños y número de 
niveles de gris sean potencias enteras de 2. Por ejemplo, un tamaño típico 
comparable en calidad a una imagen de TV monocromática es el de una 
matriz de 512 x 512 con 128 niveles de gris. 


Una de las características que define el procesado de imágenes es el ele- 
vado volumen de información que conlleva este tipo de procesos. Este he- 
cho hace que el coste computacional de las operaciones que se efectúen 
con imágenes sea muy elevado, a menos que éstas sean muy sencillas. 


—— ELEMENTOS DE UN SISTEMA DE PROCESADO 
DIGITAL DE IMAGENES 


Los elementos de un sistema de procesado de imagen básico y de pro- 
pósito general son los de la figura 1. 


| 


Almacenamiento 


NS Ordenador 
Digitalizador| — Digital => 


Imagen E Consola de 
Operación 


Visualizador 


Fig. 1. Elementos de un sistema de procesamiento digital de imágenes. 


La operación de un sistema de este tipo puede dividirse en tres catego- 
rías principales: digitalización, procesado y visualización. 


a) Digitalizadores 


Convierten una imagen en su representación numérica, de forma que 
pueda ser introducida en un ordenador digital. Los dispositivos digitaliza- 
dores más comunes son: 


— Microdensitómetros. En ellos la transparencia o fotografía se 
monta en una plataforma sobre la que se lanza un fino haz de luz. La pla- 
taforma se hace girar o mover y se van obteniendo los valores puntuales. 
El rayo reflejado (o el que la atraviesa, en el caso de transparencias) es en- 
focado hacia un fotodetector y el nivel de gris es codificado en base a la 
intensidad del rayo en cada punto. Una imagen digital se obtiene permi- 
tiendo sólo valores discretos de intensidad y posición a la salida. Son dis- 
positivos lentos pero de una gran precisión. 


— Scanners «Flying Spot». También operan sobre el principio de 
enfocar un haz de luz reflejado o transmitido sobre un foto-detector. En 
este caso, sin embargo, la imagen permanece estática y la fuente luminosa 
es un tubo de rayos catódicos (CRT) en el cual un haz de electrones incide 
sobre una superficie fluorescente. El rayo produce un punto de luz que se 
mueve y analiza la imagen sobre la superficie del tubo. Ya que el rayo se 
mueve electrónicamente, permite grandes velocidades de lectura. Tam- 
bién permite gran flexibilidad a la hora de leer una imagen, ya que la po- 
sición del rayo puede ser rápida y fácilmente establecida externamente. 


— Disectores de imagen. Su operación se basa en el principio de 
emisión electrónica. 


— Cámaras de TV. Muchos digitalizadores de cámaras de TV em- 
plean un tubo vidicón, cuya operación se basa en el principio de la foto- 
conductividad. Una imagen enfocada en la superficie del tubo produce un 
patrón de variación de conductividad que corresponde a la distribución 
de intensidad de la imagen óptica. Estos dispositivos son menos precisos 
que los anteriores, pero generalmente más baratos. 


b) Procesadores de imagen 


Los dispositivos usados para el procesado de imagen varían desde mi- 
croprocesadores para aplicaciones de propósito general hasta enormes sis- 
temas de computadores capaces de realizar gran variedad de funciones en 
imágenes matriciales de alta resolución (computadores de propósito espe- 
cial). 
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c) Visualizadores 


Son los dispositivos que se encargan de convertir las matrices numéri- 
cas almacenadas en el ordenador en una representación adecuada para la 
interpretación humana. Hay tres tipos principales de dispositivos visualiza- 
dores: 


— C.R.T. Los tubos de rayos catódicos (C.R.T.) son sistemas en los 
que la posición horizontal y vertical de cada elemento de la matriz imagen 
es convertida en voltajes que se usan para dirigir un haz de electrones, de 
forma que se produce la imagen de salida bidimensional. En cada punto 
de la pantalla, la intensidad del haz es modulada por un voltaje proporcio- 
nal al valor del punto correspondiente en la matriz (pixel). 


— Monitores de TV. Estos dispositivos convierten la imagen alma- 
cenada en el ordenador en una señal de vídeo que puede ser visualizada. 
Su ventaja es que pueden ofrecer una gama de tonalidades que se aproxi- 
ma en calidad a la de una placa fotográfica, produciendo una salida que 
puede ser fácilmente asimilada por el sistema visual. Su inconveniente es 
que el monitor debe tener un ciclo de refresco de la imagen con una fre- 
cuencia de unos 30 «frames» (marcos o imágenes) por segundo para evi- 
tar el parpadeo, mientras que la mayoría de los ordenadores convencio- 
nales no son capaces de transferir datos a esta velocidad. Para solucionar 
esto, se dota al sistema de TV de un medio de almacenamiento o «buffer» 
para transferir datos al monitor a la frecuencia de vídeo. 


— Impresoras. Los dispositivos de impresión de imagen son útiles 
para trabajos de procesado de imágenes a baja resolución. Un medio de ge- 
nerar los diferentes niveles de gris de la imagen directamente en el papel 
es aprovechando la capacidad de «over-strike» (sobre-impresión) de las im- 
presoras convencionales. El nivel de gris de cada pixel puede ser contro- 
lado por el número y la intensidad de los caracteres sobreimpresionados 
en él. Seleccionando adecuadamente el conjunto de caracteres, puede con- 
seguirse una distribución de grises aceptable con un programa sencillo y 
relativamente pocos caracteres. Este sistema de sobreimpresión de carac- 
teres es el que se ha empleado para visualizar las imágenes definidas en 
los programas del apéndice 2, proporcionando un grado de resolución su- 
ficiente para el propósito de entrenamiento al que se hallan dirigidos. 
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OPERADORES SOBRE IMAGENES — 


MUESTREO 
Y CUANTIFICACION DE IMAGENES 


ARA poder procesar una imagen f(x,y), esta función debe 
pasar por un proceso previo de digitalización tanto espa- 
cial como en amplitud. La digitalización de las coordena- 
das espaciales es lo que se llama muestreo de la imagen 
y la digitalización en amplitud (cuantificación en niveles 
de grises). 

Una imagen continua la podemos aproximar por una 
matriz de dimensiones N Xx N, de forma que: 


f(0,0) f(0,1) e f(0,N-1) 
f(x,y)= £(1,0) f(1,1) e f(1,N-1) 
f(N-1,0) f(N-1,1) a f(N-1,N-1) 


En esta matriz, cada elemento f(x,y) es un valor discreto. Una matriz 
de este tipo se denomina imagen digital, y cada uno de los elementos re- 
cibe el nombre de pixel. 

Para realizar el proceso de digitalización, hay que fijar previamente el 
valor de N, así como el número de niveles de grises que se van a permitir 
para cada pixel. Normalmente, se eligen para estos valores potencias en- 
teras de dos, de manera que: N=2" y G=2", siendo G el número de niveles 
de gris. Según esto, el número de bits necesario para almacenar una ima- 
gen digitalizada será de N x N pixels de dimensión y m bits/pixel será: 


b=NxNxm 


Así, por ejemplo, para almacenar una imagen de 128 x 128 pixels y de 
64 niveles de gris se necesitarán 98.304 bits. 

Un factor importante a tener en cuenta es la determinación de cuántos 
niveles de gris y qué valor de N se necesitan para obtener una buena aproxi- 
mación de la imagen. La resolución depende mucho de los valores que se 
denamyamnN, de forma que cuanto más altos sean éstos, mayor será la 
resolución de la imagen. Sin embargo, los requerimientos de memoria y 
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proceso aumentan considerablemente al aumentar estos parámetros. Para 
obtener una imagen de calidad análoga a la de TV se necesitan 512 x 512 
pixels y 128 niveles de grises. Como norma general, un sistema que pro- 
cese imágenes debe ser capaz de mostrar imágenes de 256 x 256 pixels y 
64 niveles de grises, como mínimo. 

El efecto que produce la reducción de N es una imagen digital, es la 
aparición de un efecto damero en la misma. Al disminuir m, el efecto que 
se produce es la aparición de falsos contornos. Estos son los efectos que 
produce la variación de N y m por separado; en 1965 Huang estudió las 
relaciones entre los dos parámetros, haciéndolos variar a un tiempo. Para 
este experimento utilizó tres imágenes con distintos niveles de detalle; a 
las conclusiones a que llegó fueron: 


1. Al aumentar N y m aumenta la calidad de imagen. 

2. En algunos casos, para un N fijo se mejora la imagen al disminuir 
m, esto es debido a que se aumenta el contraste. 

3. En imágenes con una gran cantidad de detalle, se necesitan sólo 
unos pocos niveles de grises para conseguir una buena aproximación. 


=== MUESTREO Y CUANTIFICACIÓN NO UNIFORMES 

Fijado el valor de N, en muchos casos es posible mejorar el aspecto de 
la imagen usando un método adaptativo, en el que el proceso de muestreo 
depende de las características de la imagen. En general, se necesita un ma- 
yor número de muestras en aquellas regiones en que las transiciones en- 
tre distintos niveles de grises son más pronunciadas, haciéndose, por tan- 
to, un muestreo más preciso, pudiéndose tomar un número menor de 
muestras en las regiones de la imagen donde las transiciones son más 
suaves. 

Vamos a considerar, por ejemplo, la imagen de un rostro sobre un fon- 
do uniforme. El fondo queda representado adecuadamente si tomamos un 
número pequeño de muestras, mientras que no ocurre lo mismo con el ros- 
tro, por lo que hay que aumentar el número de muestras tomadas de este 
último, mejorándose así la imagen. 

Se ve claramente que hay imágenes en las que distribuir el número de 
muestras de forma no uniforme resulta ventajoso. Una buena técnica con- 
siste en tomar un número mayor de muestras en las fronteras de transi- 
ción entre niveles de grises. 

El muestreo no uniforme se usa para resaltar los bordes; es por esto 
que no debe usarse en imágenes con regiones uniformes relativamente pe- 
queñas, por ejemplo, la imagen de una masa de gente apiñada. También 
se usa esta técnica cuando se quiere obtener una imagen más detallada. 
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Se usa una técnica parecida a la expuesta anteriormente en el mues- 
treo no uniforme, pero jugando aquí con la distribución de los niveles de 
grises. 

Hay varias formas de realizar una cuantificación no uniforme sobre una 
imagen. Una de ellas consiste en utilizar pocos niveles en los bordes, de 
manera que se consiga resaltarlos, y tomar más niveles en los cambios sua- 
ves, consiguiendo la reducción y en algunos casos la desaparición de fal- 
sos contornos. 

Otra alternativa consiste en realizar una mayor cuantificación en los ni- 
veles de grises que sean más frecuentes, haciendo una cuantificación me- 
nor en aquellos niveles que sean menos frecuentes. 

Al igual que el muestreo no uniforme, la cuantificación no uniforme se 
usa para resaltar bordes y obtener imágenes con mayor detalle de represen- 
tación. 
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LA TRANSFORMADA DE FOURIER 


Uno de los problemas que plantea el estudio de los fenómenos físicos 
es el de su representación. Dicha representación requiere la incorporación 
de conceptos que son idealizaciones introducidas por el observador como 
modelización de la realidad. Las representaciones más habituales de una 
señal son la temporal (o espacial) y la de frecuencias temporales (o espa- 
ciales). Esta última transformada de Fourier de la primera, permite acce- 
der a las componentes espectrales de la señal (representación armónica). 
Ambos tipos de representación son completas y totalmente equivalentes. 


Como es bien conocido, una función periódica puede ser descompues- 
ta mediante suma de infinitos senos y cosenos con diferentes amplitudes, 
frecuencias y fases. En el caso de funciones no periódicas, la serie infinita 
se convierte en una integral de la forma: 


+00 


F(u) = | f(x) exp(=j2rtux) dx (3,1) 


en donde exp(-j2rux) = cos(2rux) - jsen(2mux), siendo ¡=-1 la unidad 
imaginaria. 
La Ec. (3.1) proporeiona la transformada de Fourier unidimensional 


F(u) de la señal f(x). Si la función f(x) es real, su correspondiente transfor- 
mada de Fourier es compleja, es decir: 


F(u)= R(u) + jI(u) (3.2) 
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DOMINIO ESPACIAL DOMINIO DE FOURIER 


Y 
f(x,y) => 


pd Y 
£, 


(c) 


Fig. 1. Representación en el dominio espacial y de frecuencia de funciones sinu- 
soidales (a,b), bidimensionales y de una función rectángulo bidimensional (c). 
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Fig. 2. Módulo de la transformada de Fourier (b) de imágenes texturadas (a). 
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siendo R(u) e I(u) las partes real e imaginaria de F(u). También puede re- 
presentarse F(u) en la forma módulo-argumental, es decir: 


F(u) = M(u) expGD(u)) . (3.3) 
siendo M(u) = VR? + P 


y D(u) = tg”(I(u)/R(u)) 

Al módulo de la transformada de Fourier F(u) se le denomina espectro 
de Fourier y a DP(u) ángulo de fase. En la figura 1 se presentan algunos 
ejemplos de señales con variación sinusoidal y rectangular (izquierda) y 
su correspondiente transformación de Fourier (derecha). En la figura 2 se 
representan los módulos de la transformación de Fourier de algunos ejem- 
plos de imágenes con información de textura. 
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== Transformada de Fourier discreta 


Sea f(x) una función continua que se ha discretizado en la secuencia 
£(x0), K(X, + 8), f(x, + 28), ..., £(x, + [N-1]5)!, tomando N muestras. separa- 
das por Ú unidades. 

Definimos ahora f(x) = f(x,+x0), y la variable x tomará valores discre- 
tos 0,1,2,...,N-1, es decir, consideramos la secuencia (f(0), £(1), ..., K(N-1)! 
para representar los valores de N muestras equidistanciadas, tomadas de 
una función continua f(x). 


A 


Con esta notación, la transformada de Fourier discreta se define como: 


N-1 


F(u) = 1/N y f(x) exp[-321ux/N], 


x=0 
para u=0, 1,2,..., N-1, y: 
N-i 
f(x) = Y, Fu) explj21ux/N], 
u=0 


para x=0,.1, 2, «, N-1 


Los valores u =0, 1, ..., N-1 de la transformada continua se correspon- 
den con las muestras de la discreta para 0, 6, 26, ..., (N-1)6. Por otro lado, 
la relación que liga 5(x) con S(u) viene dada por la siguiente expresión: 


S(u) = 1/N S(x) 
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Extensión caso bidimensional 


Para dos variables la transformada de Fourier discreta es: 


F(u,v) = 1/MN y y f(x, y) exp[-j27(ux/M + vy/N)] 


x=0  y=0 


parau=0, 1,..,M-1yv=0,1,.., N-1 


M-1 N-1 
y f(x,y) = A by F(u,v) exp[¡2xr(ux/M + vy/N)] 


u=0  v=0 
parax=0,1,..,M-1ley=0,]1,..., N-1 
Igual que en el caso unidimensional, f(x,y) representa muestras de la 
función f(x, + x9x, y, + y0y) para x=0, 1, ..., M-1ey=0, 1,..., N-1. Aná- 
logamente también las relaciones entre los incrementos en los dominios 
del espacio y la frecuencia son: 


Su = 1/MÓ6x 
Ov = 1/NOy 


Cuando se muestrea una señal, se obtiene una matriz cuadrada; por tan- 
to, M=N, y entonces: 


N-1 


F(u,v) =1/N Y) Y f(x.y) exp[-j2r(ux + vy)/N] 
y=0 


x=0 
para u,v=0, 1,..., N-1, y: 
N-1  N-1 
f(x,y) = 1/N Y f(u,v) exp[j2xr(ux + vy)/N] 
u=0  v=0 
para x,y =0, 1,...., N-1 
Hay que observar en estas operaciones que se ha incluido el término 


1/N, puesto que es una constante multiplicativa que no influye en el par 
F(u,v) y f(x,y) transformadas. 
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Propiedad de separabilidad 


La transformada discreta de Fourier para el caso bidimensional, se pue- 
de expresar como: 


N-1 N-1 
F(u,v) = 1/N Y, exp (¡21ux/N) D, f(x,y) exp(-¡2rvy/N) 
x=0 y=0 
para u,v=0, 1, ..., N-1, y: 
N-1 N-1 
f(x,y) = 1/N A exp (¡2mux/N) y F(u,v) exp(¡2rvy/N) 
u=0 v=0 


para x,y =0, 1,...., N-1 


La ventaja de esta propiedad, es que F(u,v) o f(x,y) se pueden obtener 
en dos pasos por aplicación sucesiva de la transformada de Fourier unidi- 
mensional o su inversa. Esto se ve más claro si hacemos la siguiente sepa- 
ración de términos: 

N-1 
F(u,v) = 1/N Y, F(x,v) exp[-¡21ux/N] 


x=0 


en donde, e 
F(x,v) = N[1/N DY, f(x y) exp(¡21wy/N)] 
y=0 


En esta expresión, para cada valor de x, lo que nos encontramos entre 
corchetes es una transformada de Fourier unidimensional, con v=0, 1,..., 
N-1. Por tanto, la función bidimensional F(x,v) se obtiene haciendo una 
transformada unidimensional por filas sobre f(x,y) y multiplicando el re- 
sultado por N; por último, para obtener F(u,v), se realiza otra transforma- 
da unidimensional por columnas sobre F(x,v). 

Nótese que se obtendrían los mismos resultados haciendo primero la 
transformada por columnas sobre f(x,y) y a continuación haciéndola por 
filas sobre el resultado. 


Jl 


Representación de la transformada de Fourier 


Los espectros de muchas imágenes decrecen al aumentar la frecuen- 
cia, por lo que sus términos para frecuencias altas se hacen impercepti- 
bles. Para compensar esto, y lograr una mejor visualización, se utiliza una 
representación logarítmica del tipo: 

D(u,v) = log[1 + K abs(F(u,v))] 
en lugar de F(u,v). La función D(u,v) no es nunca negativa. 


| 


| 
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Con esto se pretende que en el espectro de la función que representa 
una imagen todos los picos tengan una amplitud similar, de forma que se 
consiga el mismo brillo e intensidad en su representación.Vamos a demos- 
trarlo viendo que, en efecto, la relación entre el máximo y el mínimo va- 
lor que puede tomar la función es menor en D(u,v) que en F(u,v). Sea: 


R == log(1 5 Emnax)/108(1 + Eo) 


donde F,,,. es el valor máximo alcanzado por F(u,v), y F.,¡, es el mínimo 
de los valores que alcanza. 

La relación Fma/Fmih es siempre mayor que R. Por ejemplo, sea 
max = 1.0, y Fmin = 0.2, tenemos entonces EFmnax/Fmin = 5.0, y R = 3.8; si ade- 
más la función F(u,v) la multiplicamos por una constante k para aumen- 
tarla, la relación F..,/Fmi, sigue valiendo 5.0, mientras que R, si k = 20, 
toma el valor 1.9, que es mucho menor. 


ls! 


IN 


Transformada rápida de Fourier 


El número de multiplicaciones y sumas complejas que requiere el 
cálculo de la transformada discreta de Fourier de la expresión: 
N1 


F(u) = 1/N Y, f(x) exp[-j210ux/N] 


x=0 


es proporcional a N?, ya que son necesarias N multiplicaciones complejas 
de la expresión f(x)exp[-¡21ux/N] y (N-1) sumas. Por otro lado, los valo- 
res de la exponencial compleja pueden calcularse una sola vez y almace- 
narse en una tabla para posteriores cálculos. 

Se ha propuesto un procedimiento de descomposición denominado 
transformada rápida de Fourier (FFT), que reduce el número de multipli- 
caciones y sumas a Nlog,N, lo cual representa un coste computacional mu- 
cho menor que el que requiere la transformación directa, en especial cuan- 
do N es grande, como puede comprobarse en la siguiente tabla: 


N N? Nlog,N 
2 4 2 
4 16 8 
8 64 24 
16 256 64 
32 1024 160 
64 4096 384 
128 16384 896 
256 65536 2048 
512 262144 4608 
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La mayor aplicación de las transformadas está en el diseño de filtros 
para mejora de imágenes y también en algunas técnicas de compresión de 
imágenes. En el apéndice 1 se describe con más detalle el algoritmo que 
permite el cálculo de la FFT. Asimismo, este algoritmo puede ser utilizado 
para evaluar la transformada de una imagen 2-D, efectuando, en primer lu- 
gar, una transformación por filas y posteriormente por columnas. En for- 
ma de pseudocódigo sería: 


ALGORITMO FFT2D 
FOR i:=1 TO N-1 DO 
BEGIN 
X(1):=£(i,j) 
FFT(N,X) 
f(i,j):=X 
END; 
FOR j:=1 TO N-1 DO 
BEGIN 
X():=fGj) 
FFT(N,X) 
f(ij):=X 
END. 


Cada llamada para calcular la transformada rápida de Fourier compor- 
ta un coste computacional proporcional a Nlog,N. Puesto que existen 2N 
de tales llamadas, el coste computacional total de la FFT bidimensional es 
proporcional a N?log,(N?). 

Este procedimiento supone un gran avance para la computación de la 
transformada de Fourier, sobre todo cuando N es un número muy grande. 
Por ejemplo, una implementación directa de la transformada de Fourier 
para N = 8192, necesita del orden de 45 minutos en una máquina como 
un IBM 7094; el mismo trabajo realizado con la implementación de la 
transformada rápida de Fourier emplearía cinco segundos en la misma má- 
quina. 
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TECNICAS DE PROCESADO 
DE IMAGENES 


AS principales técnicas de procesado de imágenes pueden 
dividirse en dos grandes categorías: métodos en el domi- 
nio de la frecuencia y métodos en el dominio espacial. Los 
primeros se basan en efectuar la transformada de Fourier 
de una imagen, mientras que los segundos manipulan di- 
rectamente los pixels de la imagen. 

El fundamento de las técnicas en el dominio de la fre- 
cuencia es el teorema de convolución. Sea g(x,y) una ima- 
=== gen formada por la convolución de una imagen f(x,y) y un 
operador invariante espacialmente h(x,y), es decir: 


g(x,y) = h(x,y) * f(x,y) 


Entonces, por el teorema de convolución, tenemos la siguiente rela- 
ción en el dominio de la frecuencia: 


G(u,v) = H(u,v) F(u,v) 


donde G, H y F son las transformadas de Fourier de g, h y f, respectiva- 
mente. La transformada H(u,v) es a veces llamada función de transferen- 
cia del proceso. Ambas ecuaciones son muy utilizadas en mejora y restau- 
ración de imágenes, siendo más importante la segunda, ya que la convo- 
lución discreta es a menudo realizada más eficientemente en el dominio 
de la frecuencia, mediante un algoritmo FFT. 

En un problema típico de mejora, f(x,y) es dada y el objetivo, tras 
calcular F(u,v), es seleccionar un H(u,v) tal que la imagen deseada, dada 
por: 


- g(%y) = F” [H(u,v) F(u,v)] 


exhiba alguna característica importante de f(x,y). Por ejemplo, los bordes 
en la imagen f(x,y) pueden ser acentuados utilizando una función H(u,v) 
que enfatice los componentes de alta frecuencia de F(u,v). 
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Las técnicas en el dominio espacial se basan en transtormaciones de 
los niveles de gris, donde el tipo de operación utilizada usada depende del 
criterio empleado para el realce o mejora, operando siempre sobre el va- 
lor de los pixels de la imagen. . 

Por otra parte, podemos hacer otra división de los métodos de mejora 
de imagen en el dominio espacial atendiendo al número de pixels de la 
imagen original, de los que depende el valor de cada pixel de la imagen 
resultante: 

— Operaciones de punto o pixel individual: En ellas el valor de cada 
pixel de la imagen resultante es función unicamente del valor que tenía 
ese mismo pixel en la imagen original, mediante una transformación de al- 
gún tipo. Dentro de este grupo se encuentran las técnicas de modificación 
del contraste por ecualización del histograma de niveles de gris. 

— Operaciones de área: En ellas se calcula el valor de cada pixel como 
una función de los valores de los pixels de un cierto entorno de dicho pixel 
en la imagen original y no sólo del valor original de ese pixel. En esta ca- 
tegoría se encuadran las operaciones de convolución y algunas transfor- 
maciones geométricas. 


0 


MANIPULACION DEL CONTRASTE 


La manipulación del contraste es una transformación radiométrica 
puntual (pixel a pixel), cuyo objetivo es lograr una mejor discriminación 
de imágenes con bajo contraste. Cada nivel de gris de la imagen se modi- 
fica mediante la transformación seleccionada, sin considerar los valores 
de los niveles de gris de los pixels adyacentes. El análisis de los histogra- 
mas de una imagen constituye en la mayoría de las ocasiones un paso pre- 
liminar para lograr una eficiente manipulación del contraste de una ima- 
gen. El histograma de una imagen proporciona el porcentaje del número 
total de pixels que presentan un valor determinado de nivel de gris. El con- 
cepto es análogo al de densidad de probabilidad que se utiliza en estadís- 
tica. Las características generales de los histogramas de niveles de gris en 
el caso de algunos tipos de imágenes aparecen en la figura 1. 


Escena con baja radiancia Escena con alta radiancia 


% % 


Nivel de gris Nivel de gris 
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Escena con bajo contraste Escena con alto contraste 


ML 


Nivel de gris Nivel de gris 


Fig. 1. Características de los histogramas de diferentes tipos de escenas. 


El histograma de una imagen sólo especifica el número total de pixels 
de cada nivel de gris y no proporciona información acerca de la distribu- 
ción espacial de los niveles de gris. La figura 2 presenta el resultado de la 


aplicación de una transformación de los niveles de gris en el histograma 
de una imagen. 


NG | | 


Fig. 2. Resultado de la aplicación de una transformación lineal en el histograma. 


Cada nivel de gris de la imagen de entrada es proyectado en la curva 
de transformación, proporcionando un nuevo valor en la imagen de sali- 
da, conservándose el área del histograma obtenido. Esta simple operación 
lineal de mejora del contraste se denomina compresión del contraste. En 
ocasiones, y con objeto de conseguir un mejor contraste en la imagen de 
salida, se suele considerar algún tipo de saturación en ambos extremos del 
rango de salida, como se muestra en la figura 3. 


Si el histograma es asimétrico, como sucede frecuentemente, es real- 
mente complicado poder controlar mediante una simple transformación 
lineal tanto el valor medio del nivel de gris de la imagen de salida así como 
la cantidad de saturación en los extremos de la transformación efectuada. 
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% 


e 
> 


t t NG 


Fig. 3. Mejora lineal del contraste con saturación variable. 


Es en estos casos en donde podría aplicarse una transformación lineal por 
tramos como la indicada en la figura 4. 


Pao a t 1 NG NG' 


Fig. 4. Mejora del contraste no lineal. 


Una operación de manipulación del contraste más general lo constitu- 
ye la denominada ecualización del histograma, que se utiliza en una muy 
numerosa cantidad de imágenes. En este caso, se utiliza como función de 
transformación la distribución acumulativa de niveles de gris (el número 
total de pixels en el histograma entre cero y cada nivel de gris). La ecua- 
lización del histograma produce una reducción del contraste en áreas muy 
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claras o muy oscuras, expandiendo los valores intermedios. Un ejemplo de 
ecualización del histograma está representado en la figura 5. 


1 ! NG” 


Fig. 5. Ecualización del histograma. 


Entre los tipos de operaciones que pueden incluirse dentro de las de 
manipulación del contraste se encuentra la denominada como umbraliza- 
ción. El objetivo de la umbralización no es mejorar el contraste de la ima- 
gen considerada, sino segmentarla en dos clases diferenciadas: objetos y 
fondo, mediante la simple aplicación de un umbral en el rango de los ni- 
veles de gris. Para llevar a cabo la operación de umbralización suponga- 
mos que los niveles de gris de una imagen dada f(x,y) presentan el histo- 
grama de la figura 6. 


B1 : B2 
NEGRO T ; BLANCO 


Fig. 6. Umbralización del histograma. 


De la observación de la figura puede verse que la mayoría de los pixels 
de f(x,y) se encuentran concentrados en la zona más oscura, estando el res- 
to repartidos más o menos uniformemente. Este histograma es caracterís- 
tico de las imágenes que contienen objetos grises en un fondo oscuro. Para 
hallar la frontera entre los objetos y el fondo, dividiremos el histograma 
en dos zonas separadas por un umbral T. Se trata de seleccionar T de tal 
forma que la banda B1 contenga los niveles asociados al fondo, mientras 
que la banda B2 contenga los niveles de los objetos. Así, un cambio en el 
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nivel de gris de una banda a otra indicará la presencia de una frontera. 
Para conseguir la detección de fronteras tanto en sentido horizontal como 
vertical, será necesario explorar la imagen por filas y por columnas. Así 
pues, los pasos a seguir serían los siguientes: , 


a) Para cada fila de f(x,y), crear una fila de una imagen intermedia 
g1(x y), utilizando la siguiente relación para y=0,1,...,N-1: 


LE si los niveles de gris de f(x,y) y de f(x,y-1) g1(x,y) están en diferentes 
bandas 
LB si se encuentran en la misma banda 


b) Para cada columna de f(x,y), crear una columna de una imagen in- 
termedia g2(x,y), utilizando la siguiente relación para x=0),...,N-1: 


LE si g1(x,y) o g2(x,y) son iguales a LE g2(x,y) 
LB si no son iguales a LE 


SUAVIZADO DE IMAGEN 


Las operaciones de suavizado se usan principalmente para reducir efec- 
tos espúreos que pueden presentarse en una imagen digital como conse- 
cuencia de un sistema de muestreo pobre o del canal de transmisión. Pue- 
den realizarse tanto en el dominio espacial como en el de Fourier. 


a) En el dominio espacial. Promedio de un entorno. 


Dada una imagen f(x,y) se intenta generar una imagen suavizada g(x,y) 
cuyos niveles de gris en cada punto (x,y) se obtienen promediando los ni- 
veles de gris de los pixels de f contenidos en un entorno predefinido de 
(x,y). En otras palabras, la imagen suavizada se obtiene usando la relación: 


8(x,y) = ar E f,m) V ¡(n,m) eS 


para x,y = 0,1,2,...,.N-1. S es el conjunto de coordenadas de vecinos del pun- 
to (x,y), pero sin incluir a éste, y M es el número total de estos puntos. 

El tamaño del entorno elegido es arbitrario, pero hay que tener en 
cuenta que el proceso de suavizado de bordes puede producir un efecto 
de «blurring» o emborronamiento en la imagen, que no es deseado. El gra- 
do de emborronamiento será, además, fuertemente proporcional al radio 
del entorno usado. Para un radio dado, el efecto de emborronamiento pue- 
de ser reducido usando un procedimiento de umbralización. Es decir, en 
lugar de usar la ecuación anterior para calcular g(x,y) en todos los puntos, 
usaremos dicha expresión si 


f(x, y) — a 2 mn) > T 
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y en caso contrario diremos que g(x,y) = f(x,y); donde T es un umbral no 
negativo predefinido. De esta forma podemos reducir el ruido de la ima- 
gen evitando en gran parte el efecto indeseado de «blurring». 


b) En el dominio de Fourier. Filtros paso bajo. 


Los bordes y otras transiciones pronunciadas (como el ruido) en los ni- 
veles de gris de una imagen, contribuyen grandemente al contenido de alta 
frecuencia de su transformada de Fourier. Por tanto, el emborronamiento 
podría ser conseguido, vía dominio de la frecuencia, atenuando un deter- 
minado rango de componentes de alta frecuencia en la transformada de 
una imagen dada. 


Recordemos la ecuación: 
G(u,v) = H(u,v) F(u,v) 


donde F(u,v) es la transformada de la imagen que queremos suavizar. El 
problema es determinar una función H(u,v) que nos produzca G(u,v) por 
atenuación de las componentes de alta frecuencia de F(u,v). La transfor- 
mada inversa de G(u,v) nos dará posteriormente la imagen suavizada de- 
seada g(x,y). Ya que las componentes de alta frecuencia son «eliminadas» 
O filtradas y la información en el rango de bajas frecuencias es mantenida 
sin atenuación, este método se conoce comúnmente como filtrado de pa- 
so-bajo. La función H(u,v) se denomina función de transferencia del filtro 
considerado. En todos los casos que se considerarán más adelante, las ope- 
raciones de filtrado afectan de la misma manera a la parte real y a la parte 
imaginaria de la transformada de Fourier. 
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Ejemplos de filtrados paso-bajo 


Los ejemplos siguientes corresponden a filtros paso-bajo con simetría 
rotacional, es decir, el filtro completo se generaría rotando 360” la sec- 
ción mostrada en cada uno de los ejemplos señalados. 


a) Rectangular. Un filtro paso-bajo cuya función de transferencia sa- 
tisfaga la relación 


_[1siR=<RO 
a CE 


Gráficamente, en la figura 7 se representa el perfil de este tipo de filtro, 
que suele denominarse «ideal», puesto que deja pasar todas aquellas fre- 
cuencias situadas dentro de un círculo de radio menor que R, sin atenuar- 
las, eliminando todas aquellas frecuencias situadas en el exterior del 
círculo. Al valor R, se le denomina frecuencia de corte, que constituye un 
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H(R) 


Ro R 


Fig. 7. Filtro paso-bajo ideal. 


concepto muy útil a la hora de definir las características de un filtro, per- 
mitiendo la comparación con otros tipos de filtros. 

El resultado de la aplicación de este tipo de filtros produce un emborro- 
namiento de la imagen de partida como consecuencia de la eliminación 
de las altas frecuencias, o, lo que es lo mismo, una pérdida de definición 
de los bordes de la imagen. 


b) Filtro de Butterworth 


La función de transferencia de un filtro de Butterworth de orden n y 
de frecuencia de corte R,, viene dada por la siguiente expresión: 


1 
AUR TA 


Gráficamente, en la figura 8 se presenta el perfil de este filtro, pudiendo 
observarse cómo esta función de transferencia no presenta un cambio tan 
abrupto como en el anterior caso de filtro ideal, siendo en general más 


H(R) 


0 1 2 3 R 


Fig. 8. Filtro paso-bajo de Butterworth. 
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aconsejable la utilización de este tipo de filtros que presentan una caída 
más suave, puesto que, como puede demostrarse, la utilización del caso 
ideal produce falsos contornos en la imagen filtrada. Seguidamente vere- 
mos algunos otros ejemplos de filtros que presentan una caída suave de la 
respuesta en frecuencias. 


c) Filtro trapezoidal 


Este tipo de filtro constituye un compromiso entre el filtro ideal y un 
filtro con caída suave. La función de transferencia de este filtro viene dada 
por la siguiente expresión: 


1 si R < RO 


1 : 
HR) = | RORTÓRRDO siRO<R=<RI 
0 si R= RI 


Gráficamente, el perfil de este tipo de filtro viene representado en la figu- 
ra 9. 


H(R) 


0 R, R, R 


Fig. 9. Filtro paso-bajo trapezoidal. 


d) Filtro exponencial 


Este tipo de filtro paso-bajo, que también se utiliza frecuentemente en 
procesamiento de imágenes, presenta una frecuencia de corte situada a dis- 
tancia RO y una función de transferencia dada por la siguiente expresión: 


H(R)= exp[(R/R0)”] 
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Gráficamente, el perfil de este tipo de filtros viene representado en la 
figura 10. 


H(R) : 


' 
0 1 2 3 R 


Fig. 10. Filtro paso-bajo exponencial. 


A 


REALCE DE BORDES 


Puesto que, como ya se ha señalado, los cambios abruptos de los nive- 


0 


Ejemplos de filtrados paso-alto 
a) Filtro rectangular 


Este tipo de filtros presenta una función de transferencia del tipo: 


_JOSiR<RO 
A E 


siendo R, la frecuencia de corte. Gráficamente, este tipo de filtros viene 
representado en la figura 11. 


b) Filtro de Butterworth 


La función de transferencia de un filtro paso-alto de Butterworth de or- 
den n y frecuencia de corte RO viene dada por la siguiente expresión: 


1 
H(R) = 7+TR0/R] 2n 


W 


O 


H(R) 


0 R, R 


Fig. 11. Filtro paso-alto ideal. 


Gráficamente, este tipo de filtros viene representado mediante la figu- 
ra: 12. 


Fig. 12. Filtro paso-alto de Butteworth. 


c) Filtro exponencial 


La función de transferencia de un filtro paso-alto exponencial con fre- 
cuencia de corte RO, viene dada por la siguiente expresión: 


. H(R)= exp[-RO/R]" 


en donde el parámetro n controla la velocidad de subida de la función ex- 


ponencial. Una gráfica del perfil de este tipo de filtros viene dada en la fi- 
gura 13. 
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Fig. 13. Filtro paso-alto exponencial. 


d) Filtro trapezoidal 


La función de transferencia de un filtro paso-alto trapezoidal viene dada 
por la siguiente expresión: 


0 siR<RIl 
A 1 : 
si R > RO 


en donde RO y R1 son parámetros conocidos y se cumple que RO>0R1. 
Gráficamente, este tipo de filtros viene representado en la figura 14. 


H(R) 


0 R R, R 


Fig. 14. Filtro paso-alto trapezoidal. 
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TRANSFORMACIONES GEOMETRICAS 


En ocasiones es necesario efectuar cambios geométricos en una ima- 
gen digital, como, por ejemplo, cambiar de sistema de coordenadas, efec- 
tuar rotaciones, etc. Una vez realizada la transformación de coordenadas 
es necesario efectuar un proceso de interpolación con objeto de ajustar la 
imagen resultante a la retícula x,y, en donde se encuentra situada la ima- 
gen original. Por ejemplo, la transformación de coordenadas de un cua- 
drilátero definido por cuatro puntos puede efectuarse a través de los si- 
guientes polinomios: 


x=4,+9,X +38 Y +ary +2 Ry +a yy 

y=b, +b Xx + by +bay +b (2 +b, (y) 
en donde las coordenadas con prima son las de la imagen de salida y las 
coordenadas sin prima corresponden a la imagen original. Si los coeficien- 
tes de (x”)? e (y”) son cero, los ocho restantes coeficientes pueden calcu- 
larse a partir de los cuatro puntos correspondientes a los vértices del cua- 
drilátero de entrada. Los tipos de transformaciones posibles en este últi- 
mo caso vienen representadas en la figura 15. 


DESPLAZAMIENTO ESCALADO 

x=apg+x' pre 

y =bp+y y =09 
ESC AMIENTO PERSPECTIVA 
x=x' +07 x=axy 
ye. yy 


GIRO 

x=a¡x +0) 
y =bjx' + bay” 
a¡= y cos Y 
a, =-b, =sen Y 


n 


Fig. 15. Ejemplos de transformadas geométricas. 
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VISION ARTIFICIAL 4 


== A mayor parte del desarrollo en visión artificial se ha lle- 
vado a cabo en aplicaciones industriales en donde el en- 
torno en que se encuentra la máquina es bien conocido y 
lá tarea encomendada al sistema de visión está claramen- 
te establecida. Un ejemplo típico lo constituye el sistema de 
visión utilizado por un robot para manipular objetos en 
movimiento situados en una cinta transportadora. En 
otras áreas de aplicación en donde la información a ma- 
o mipular está peor definida, el desarrollo llevado a cabo es 
comparativamente menor, como, por ejemplo, sucede en el caso del aná- 
lisis e interpretación de imágenes de rayos X. 

El objetivo de un sistema de visión artificial es analizar las imágenes y 
producir descripciones de dichas imágenes. El sistema de visión puede ser 
considerado como un elemento de un bucle de realimentación que está re- 
lacionado con la captación y análisis de imágenes, mientras que otros ele- 
mentos del sistema están dedicados a la toma e implementación de dichas 
decisiones (figura 1). 


Fig. 1. Partes constituyentes de un sistema de visión artificial, 
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Existen diferentes niveles de descripción posibles. En general, se con- 
siderará suficiente aquel nivel de descripción que permita llevar a cabo la 
tarea deseada. 


Entre las imágenes más sencillas que puede tratar un sistema de visión 
se encuentran las denominadas imágenes binarias. Una gran parte de las 
imágenes industriales que se procesan son binarias, puesto que en muchas 
ocasiones las descripciones que pueden extraerse de este tipo de imáge- 
nes son suficientes para este tipo de aplicaciones. 


Una imagen binaria viene definida por una función característica b(i,j), 
que es cero para aquellos valores correspondientes al fondo de la imagen 
y uno para los puntos correspondientes al objeto. Un ejemplo de la ima- 
gen binaria está representada en la figura 2. 


Fig. 2. Ejemplo de imagen binaria l, indicando los centros de gravedad de cada ob- 
jeto. 


Esta imagen binaria puede proceder de una operación de umbraliza- 
ción como las anteriormente descritas. La gran ventaja de las imágenes bi- 
narias es que son fáciles de digitalizar, almacenar y transmitir, puesto que 
con muy pocos bits es posible representar este tipo de imágenes. 


El procesamiento de imágenes binarias requiere tan sólo el cálculo de 
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algunas propiedades geométricas y topológicas. Veamos seguidamente al- 
gunas de ellas: 


N-1 N-1 


AREA: A= 2 2 bj) (4.1) 
i=0 j= 
N-1 N-1 
x,=1/A Y, Y» ib) (4.2) 
i=0 j=0 
POSICION: 
N-1 N-1 
y=1/A Y jbG) (4.3) 
i=0 j=0 


Las expresiones de la coordenada de la posición corresponde a lo que 
se denomina momentos de primer orden y coincide con el concepto de 
centro de gravedad de un objeto de la Mecánica. 


Así como para el cálculo de la posición es necesario utilizar un mo- 
mento de primer orden, para obtener la orientación es necesario el cálcu- 
lo de momentos de segundo orden. Supongamos que se desea calcular la 
orientación del objeto que se encuentra situado en la figura 3. Para calcu- 


Fig. 3. Orientación de un objeto binario. 
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lar el eje según el cual se encuentra orientada la figura, hay que tener en 
cuenta que este eje coincide con el denominado eje de mínima inercia que 
se utiliza asimismo en Mecánica. Por ello, habrá que minimizar la siguien- 
te expresión: 


N-1 N-1 


E= Y Y rot) (4,4) 


i=0 j¡=0 


en donde r es la perpendicular trazada desde (i,j) a la línea anterior (figura 
4). Puesto que el eje de mínima inercia pasa por el centro de grávedad, es 


o 


Fig. 4. Distancia de un punto (ij) a una lína: isen8 — ¡cos0 + O, siendo 8 la inclina- 
ción respecto al eje i y O es la perpendicular a la línea desde el origen. 


conveniente efectuar un cambio del sistema de coordenadas a uno nuevo 
que pase por dicho punto y sea paralelo al antiguo: 


Y=ij= (4.5) 
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Sustituyendo en la expresión (4.4) y teniendo en cuenta la distancia de un 
punto (i,j) a una recta definida por los parámetros 8 y 0, de la Fig. 4 se ob- 
tiene la siguiente expresión: 


N-1 N-1 


E= Y Y (i' sen 0- j' cos 0)? b(i,j) (4.6) 


i=0 ¡=0 


y desarrollando la expresión anterior, queda: 


E=(a sen20-2b sen8cos0+c cos?0 (4.7) 
en donde: 


N-1 N-1 


a= Y Y bj) (1) (4.8) 


i=0 j=0 


N-1 N-1 


b= Y, Y, dG.) (15) (4.9) 


i=0 ¿=0 


N-1 N-1 


e= 2, 2, dúj G7 (4.10) 


i=0 j=0 


Para hallar la orientación 0, del eje de mínima inercia habrá que derivar 
e igualar a cero la expresión (4.7), obteniendo como resultado: 


8,= 1/2 tang"(2b/a-c) (4.11) 


En la figura 5 se presenta un procedimiento en lenguaje LISP denomi- 
nado CENTRO, que determina el centro de gravedad de una imagen rec- 
tangular de dimensión N filas x M columnas, devolviendo los valores 1, y 
J, de las coordenadas del centro de gravedad. En la figura 6 se presenta 
un procedimiento LISP que devuelve una lista que contiene el centro de 
gravedad, así como el valor del ángulo que forma el eje de mínima inercia 
con el eje (1). 


sión N x M. 


O (CAR C-0-A))> (JO (CADR C-0-A)) 
(LET ((10 (CAR C-O 


Fig. 6. Función LISP que permite el cálculo de la orientación de una imagen bina- 
ria N x M. 
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Otras características geométricas de interés pueden ser el perímetro de 
la figura estudiada, que puede estimarse como el número de ceros que se 
encuentran contiguos con unos. En la figura 7 se presenta un procedimien- 
to que permite calcular el perímetro de una imagen binaria, si bien cabe 
señalar que el valor obtenido es una sobreestimación del valor real, debi- 
do a que el proceso de cuantificación produce en la frontera de los obje- 
tos un efecto en diente de sierra. 


(M_ ARRAV=-DIMENSTON IMAGE 4)) 
- ENDION 1MADE ] 


(GETO SUM (+ Slim 1)))))))) 
JET ION 1+ SU 172722722) 


dimensión N x M. 
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SISTEMAS DE PROCESAMIENTO 
- DE IMAGENES 


CARACTERISTICAS GENERALES 


OS sistemas de procesamiento de imágenes se definen como 
aquellos capaces de realizar alguna de las siguientes opera- 
ciones: 


a) Adquirir imágenes de uno o más sensores, represen- 
tándose éstas en forma de imagen. La representación de una 
imagen se efectúa mediante una matriz bidimensional, cons- 
tituyendo el nivel de gris de cada pixel la tercera dimensión. 

b) Almacenar imágenes de manera que puedan ser pro- 
cesadas. 

c) Realizar funciones de procesado tales como mejora, restauración, 
codificación, análisis, etc. 


Los componentes necesarios que requieren los sistemas de procesa- 
miento para llevar a cabo sus funciones son los siguientes: 


— SENSORES, como, por ejemplo, cámaras de TV vidicón, cámaras 
de TV de estado sólido CCD, etc. 


— INTERFACES DE CAPTACION DE IMAGENES (FRAME GRAB- 
BERS). 


— MEMORIAS DE IMAGEN, para almacenamiento y presentación. 


— PROCESADORES DE IMAGEN, incluyendo los denominados proce- 
sadores puntuales y los procesadores espaciales (array processors). 


— MONITORES. 
— MEMORIAS DE ALMACENAMIENTO. 
— INTERFACES DEL USUARIO CON EL SISTEMA. 


— SOFTWARE BASICO, que permita al usuario la realización de ope- 
raciones de procesado. 
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Descripción de los componentes del sistema 


x/ 


Sensores 


La tecnología estándar para la imagen electrónica es la cámara vidi- 
cón. La mayoría de las cámaras vidicón incluyen la electrónica necesaria 
para producir una imagen de TV convencional. La sensibilidad espectral 
del vídeo va desde 350 a 700 nanómetros. La luz es enfocada en un obje- 
tivo fotosensible que genera un modelo de imagen como carga, objetivo 
que es explorado por un haz electrónico produciendo la señal de vídeo. El 
haz electrónico se genera en un tubo de vacío, lo cual hace que este tipo 
de dispositivos sean pesados y frágiles, requiriendo tensiones elevadas 
(900 V). 

En las cámaras que utilizan tecnología de estado sólido (CCD, siglas de 
Charged Coupled Devices), la luz que incide sobre cada elemento genera 
una carga proporcional a la intensidad lumínica, recogiendo esta carga 
mediante condensadores, y transferida en forma de paquetes de carga a 
un amplificador. La imagen aparece entonces como una serie de voltajes 
a la salida del amplificador. Las ventajas más importantes de las cámaras 
CCD son su gran robusted mecánica, su alta sensibilidad, difícil saturación, 
así como un peso y un tamaño reducidos. Sus desventajas más importan- 
tes son su coste elevado, así como su baja resolución. En la actualidad y 
debido a sus características, se tiende a utilizar en entornos industriales 
cada vez en mayor medida las cámaras CCD, puesto que en este tipo de apli- 
caciones la baja resolución no impide llevar a cabo las tareas deseadas. 


b) Interfaces de captación de imágenes (frame grabbers) 


En la mayoría de las ocasiones la señal de vídeo procedente de las cá- 
maras de TV es digitalizada y almacenada en el computador como una ma- 
triz bidimensional de niveles de gris. La función de un «frame grabber» es 
la de convertir la señal de entrada en un formato que pueda ser almace- 
nado y procesado. Este proceso requiere la digitalización de la señal y su 
transferencia a la memoria de imagen. Los componentes de un «frame 
grabber» son los siguientes: 


— Adquisición de la señal de entrada. 
— Sincronización con el formato de la señal de entrada. 
— Convertidor analógico-digital de alta velocidad. 


Es importante tener en cuenta el formato de las señales de entrada a 
la hora de efectuar su captación, ya que algunas señales, como las de TV, 
incorporan los sincronismos en la señal de vídeo, mientras que en otro 
tipo de señales los sincronismos aparecen exteriormente. 
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c) Memorias de imagen 


Las imágenes procedentes del sensor se almacenan en la memoria de 
imagen en posiciones consecutivas, representando cada valor almacenado 
el nivel de gris de la imagen considerada. Debido al gran volumen de in- 
formación necesario para almacenar imágenes, las memorias de imagen 
son de gran capacidad. Así, por ejemplo, una imagen de 512 x 512 pixels 
x 8 bits de cuantificación requiere 0.25 Mbytes de almacenamiento. En 
cuanto a los tiempos de acceso, es decir, el tiempo requerido para leer o 
escribir en la memoria de imagen, no debe exceder el requerido para la 
captación de la señal de vídeo. Este hecho supone que en algunos siste- 
mas de procesado de imágenes la propia memoria de refresco se utilice 
como memoria de imagen. 

La memoria de refresco constituye la parte más importante del sistema 
de visualización, ya que su diseño es crítico. En los visualizadores recien- 
tes, el hecho de que la memoria de refresco se utilice como memoria de 
imagen, permite realizar funciones tales como ampliación, desplazamien- 
to, corte, etc. Estas operaciones no destruyen los datos contenidos en la 
memoria. En la mayoría de los sistemas, la memoria de refresco se en- 
cuentra diseñada en forma de canales. Un canal típico contiene una ima- 
gen de 512 x 512 x 8 bits, lo cual equivale a refrescar una imagen mono- 
croma en pantalla. Si se utilizan tres canales en paralelo, podrían repre- 
sentarse imágenes en color. Los sistemas actuales proporcionan entre 1 y 
64 canales, existiendo asimismo sistemas que proporcionan imágenes de 
1024 x 1024 pixels y 16 bits de precisión. 

El computador central accede a la memoria de refresco mediante al- 
guno de los siguientes métodos: 1) a través de una interface que permite 
acceder al usuario de forma aleatoria por medio de una transferencia DMA, 
o 2) considerando una subregión rectangular del canal dentro del espacio 
de direccionamiento del computador, por lo que se accede a la misma 
como si se tratase de la memoria principal del computador. 


Desde mediados de los setenta, las memorias de refresco se han cons- 
truido utilizando memorias RAM. Las sucesivas mejoras en la tecnología 
RAM ha permitido incrementar la capacidad de los chips, por ejemplo, 4 
Kbits, 16 Kbits, 64 Kbits, etc. Sin embargo, existe una relación inversa en- 
tre el tiempo de acceso a este tipo de memorias y el consumo de potencia, 
porcentaje de errores y coste. Es por ello que las memorias de refresco se 
han construido a partir de dispositivos relativamente lentos. Un valor típi- 
co del tiempo de ciclo de este tipo de memorias es el de 400 nseg., que es 
mucho más lento que el requerido para el refresco de los visualizadores. 
Para obtener las velocidades de refresco necesarias, las memorias de los 
sistemas utilizan multiplexación de pixels, típicamente con un factor de 8 
en los visualizadores que utilizan un formato de 512 x 512. En general, la 
mayoría de los sistemas multiplexan cada bit de intensidad en un chip di- 
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ferente; por ello, en un diseño típico se accede a 64 chips de RAM durante 
un simple ciclo de memoria. 

La mayoría de las limitaciones de las memorias de refresco proceden 
de la utilización de la multiplexación. Este hecho supone que algunos vi- 
sualizadores no dispongan de todos los valores enteros de ampliación 
(zoom), siendo los valores 2:1, 4:1 y 8:1 mucho más fáciles de implemen- 
tar. La multiplexación también limita el acceso por columnas, puesto que 
los ocho pixels multiplexados son direccionados por filas. La consecuen- 
cia de ello es que el acceso por columnas será ocho veces más lento que 
por filas. 


d) Monitores de presentación de imágenes 


El monitor de presentación de imágenes constituye un monitor indus- 
trial de alta calidad. Generalmente, se encuentra calibrado para presentar 
una trama cuadrada de datos con una relación de aspecto 1:1. Los tama- 
ños de presentación estándar son 512 x 512 o 1024 x 1024 pixels, si bien 
también suelen encontrarse otros tamaños tales como 640 x 480, 
1024 x 1280. La señal de TV comercial utiliza velocidades de 30 (o 25) tra- 
mas/seg., con una relación de aspecto 2:1 y entrelazamiento. Sin embar- 
go, se vienen utilizando recientemente visualizadores con un refresco de 
60 (o 50) tramas/seg. y sin entrelazamiento. Este último tipo de visualiza- 
dores produce una presentación carente de parpadeos, lo cual minimiza 
la fatiga del observador. 

Por otro lado, es bastante frecuente la utilización por parte de algunos 
fabricantes de terminología ambigua, fundamentalmente en el caso del tér- 
mino resolución. La única medida de resolución aceptable consiste en co- 
locar un test de prueba en el monitor y medir ópticamente la resolución 
con un fotómetro. 

Las principales restricciones de los visualizadores se deben a la diver- 
gencia del haz de electrones y a la dispersión de los fotones. En la figura 
1 se presenta la arquitectura típica de un visualizador de un sistema de pro- 
cesamiento de imágenes. 


e) Procesadores de imagen 


Existen un gran número de arquitecturas para procesamiento de imá- 
genes, incluyendo las arquitecturas especializadas. En general, cabe seña- 
lar que los procesadores serie convencionales tipo von Neuman no son los 
más adecuados para el procesado de imágenes, sino más bien arquitectu- 
ras de tipo paralelo. 

Las principales arquitecturas para procesamiento de imágenes son las 
siguientes: 


— Procesadores serie convencionales. 
— Procesadores con paralelismo lineal (pipeline). 
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Computador Central 


Vídeo analógico 


Fig. 1. Arquitectura de un visualizador de un sistema de procesamiento de imáge- 
nes. 


— Procesadores paralelo (array processors). 
— Procesadores de propósito especial. 
— Procesadores VLSI. 


Como anteriormente se indicó, los procesadores serie no son los más 
adecuados para el procesado de imágenes. Sin embargo, presentan venta- 


jas sustanciales con el resto de procesadores, como son su bajo coste y su 
flexibilidad. 


== Métodos de iluminación 


Un aspecto importante a tener en cuenta a la hora de efectuar la cap- 
tación de una imagen por parte de un sistema de visión lo constituye el mé- 
todo empleado para la captación de dicha imagen. Cuando lo que se desea 
es captar una escena industrial es necesario, en primer lugar, decidir el nú- 
mero y la configuración de las cámaras a emplear, puesto que si se desea 
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obtener información tridimensional, será necesario el concurso de varias 
cámaras. Otro de los aspectos a tener en cuenta es el método de ilumina- 
ción empleado, puesto que una iluminación inadecuada puede producir 
sombras, brillos u otros efectos indeseados. Losrmétodos más usuales de 
iluminación de escenas son los siguientes: 


a) Iluminación difusa. Es un tipo de iluminación indirecta que condu- 
ce a una buena eliminación de sombras, permitiendo una modelización tri- 
dimensional del objeto considerado utilizando una sola cámara (figura 2). 


CAMARA 


OBJETO 


Fig. 2. Ejemplo de disposición de una iluminación difusa. 


b) Iluminación posterior. Constituye el sistema de iluminación más uti- 
lizado. Permite el reconocimiento de piezas pequeñas y ligeras, por lo que 
es un método muy indicado para aplicaciones de control de calidad y cla- 


sificación de piezas (figura 3). 
la CAMARA 
1] 
1] 


Fig. 3. Disposición de iluminación posterior. 
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c) Iluminación modulada. Este tipo de iluminación permite la extrac- 
ción de información de curvaturas, así como de su uniformidad superfi- 
cial por medio de la deformación que experimentan unas líneas proyecta- 
das sobre el objeto, mediante la utilización de una máscara conocida (figu- 
ra 4). 

FOCO DE LUZ 


CAMARA E 
REJILLA 
mm 
OBJETO 


Fig. 4. Disposición de iluminación modulada. 


d) Iluminación direccional. Este tipo de iluminación se utiliza funda- 
mentalmente para la observación de las rugosidades superficiales de un ob- 
jeto, así como para la observación de pequeñas áreas del objeto conside- 
rado, mediante la utilización de un haz coherente (láser) o no (figura 5). 


¡$ CAMARA 
17 


Fig. 5. Disposición de iluminación direccional. 
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CARACTERISTICAS HARDWARE Y SOFTWARE 
DE ESTOS SISTEMAS 


En los últimos diez años, los sistemas para procesamiento de imágenes 
han experimentado un considerable desarrollo, de tal manera que opera- 
ciones de procesado que en los primeros sistemas suponían un elevado 
tiempo de cálculo, en la actualidad se realizan prácticamente en tiempo 
real. A partir de los años setenta, estos sistemas se han venido utilizando 
fundamentalmente en investigación. La razón fundamental del interés por 
el diseño de este tipo de sistemas radica en el hecho de que la tecnología 
serie de los computadores existentes no es la más adecuada para el proce- 
samiento de este tipo de información. 

Los primeros intentos de construcción de este tipo de sistemas se efec- 
tuaron en la década de los cincuenta. Estos sistemas se caracterizaban por 
su gran tamaño, memoria limitada y mala calidad de imagen, debida al rui- 
do. El procesamiento digital de imágenes comenzó realmente a fines de 
los sesenta, con el desarrollo de computadores y minicomputadores más 
baratos que los anteriores. Un ejemplo del software de procesamiento de 
imágenes desarrollado en esa época es el VICAR, construido por el Jet Pro- 
pulsion Laboratory. Este tipo de sistemas estaban basados en un ordena- 
dor central y las tareas se atendían mediante una cola de trabajo y eran de 
un coste excesivamente caro como para poder ser comercializados. Algo 
más tarde, el Laboratorio de Procesamiento de Imágenes, de la Universi- 
dad del Sur de California, efectuaba los primeros procesamientos locales 
con ayuda de un microcomputador. 

La arquitectura propuesta a principios de los setenta, basada en un mi- 
nicomputador junto a una serie de periféricos, se ha mantenido hasta la ac- 
tualidad. Los periféricos utilizados permiten la adquisición y digitalización 
de las imágenes, su presentación y registro. Por otro lado, y con objeto de 
aumentar la velocidad de procesamiento, se les suele conectar a un com- 
putador central o incluso añadirles un procesador de array. 

Los primeros visualizadores comerciales fueron el Aerojet (posterior- 
mente denominado Comtal) y permitían la realización de manipulaciones 
radiométricas, superposición de gráficos y acceso directo a los pixels de 
la imagen mediante cursor. Estos visualizadores estaban limitados a la pre- 
sentación de las imágenes y a manipulaciones globales de la intensidad o 
del color. 

En la mitad de los setenta, se desarrollaron los primeros sistemas de 
procesamiento realmente interactivos con la incorporación de software de 
tiempo real. Así, por ejemplo, cabe citar el sistema VICAR, que posterior- 
mente dio lugar al Mini-VICAR, el sistema GE 100, de General Electric 
para procesamiento interactivo de imágenes LANDSAT, o el sistema PE- 
COS, de Electromagnetic Systems Laboratory. La evolución de las nuevas 
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tecnologías de computadores fue acompañada de una progresiva evolu- 
ción de los sistemas de presentación. Así, la incorporación de memorias 
de acceso aleatorio RAM como memorias de refresco permitió la realiza- 
ción de operaciones más complejas, como, por ejemplo, la suma y resta 
efectuadas a velocidad de vídeo; el cálculo de histogramas por hardware, 
particiones y ampliaciones de la imagen, etc. 

A mediados de los setenta, International Imaging Systems introdujo el 
denominado Modelo 70, que incorporaba un bucle de realimentación a ve- 
locidad de vídeo, que con el paso de los años se convirtieron en una ca- 
racterística común a todos los equipos, incorporando unidades aritméti- 
co-lógicas para aumentar la velocidad (como, por ejemplo, los sistemas 
Gould DeAnza, Vicom, Ikonas, etc). El procesamiento en forma de bucle 
de realimentación refleja una de las características básicas de los algorit- 
mos de procesamiento de imágenes: la realización de operaciones simples 
con imágenes puede ser aplicada iterativamente a cada uno de los pixels. 

Los visualizadores desarrollados a fines de los setenta y principios de 
los ochenta incorporaban las siguientes características: 1) operaciones op- 
cionales para la realización de operaciones específicas, tales como opera- 
ciones de convolución, y 2) arquitecturas paralelas y bucles de realimen- 
tación para la realización de algoritmos iterativos. 


Sistemas comerciales de procesamiento 
digital de imágenes 


Los sistemas comerciales de procesamiento digital de imágenes pue- 
den ser de tres tipos: 


e Sistemas de propósito general. Están constituidos por un minicompu- 
tador, junto con los periféricos necesarios para efectuar las operaciones 
de procesamiento. En otras ocasiones se comercializa sólo el procesador 
de imágenes para ser conectado a algún computador central. Ejemplos de 
este tipo de sistemas lo constituyen el sistema VICOM y el GOULD DeAnza. 
En la figura 6 aparece un esquema de la arquitectura del sistema de pro- 
cesamiento de imágenes VICOM. Como caso particular de este tipo de sis- 
temas se encuentran los denominados de propósito específico. Constitu- 
yen sistemas cerrados utilizados generalmente para aplicaciones concre- 
tas. Así, por ejemplo, cabe citar el sistema InterScan, de Control Automa- 
tion. En la figura 7 aparece un cuadro comparativo de algunos de los sis- 
temas comerciales existentes, tanto de propósito general como de propó- 
sito específico. 

e Estaciones de ingeniería (workstations). Este tipo de sistemas propor- 
cionan una potencia local de cálculo elevada, de tal manera que el usua- 
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Fig. 6. 
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Sistema de procesamiento de imagen VICOM. 
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PROCESADO VISUALIZACION, RESOLUCION 


CARACTERISTICAS 


VELOCIDAD B/N COLOR | GRISES | PIXELS 


Tiempo real, convolución espacial 
Programas en ASSEMBLER, FORTRAN, PASCAL 
Aplic. en convolución, inspec. y análisis 


Programas en BASIC y FORTRAN 
Inspección y microscopía 


Programas en BASIC y FORTRAN 


6,5 seg. . : 
Inspección y microscopía 


Tiempo real 
Programas en FORTRAN y PASCAL 
Inspección de objetos móviles, control de medid. 


Programa en FORTH 
Aplicaciones de inspección 


Programa en FORTH 
Inspección 


2 
x 
E 


Tiempo real 
Robótica, inspección 


ES 
x 
3 


Posicionamiento de componentes en inspección 
Robótica 


45 seg. Inspección visual 


Accionamiento de formas, inspección 
Robótica 


Programa en BASIC 


10 seg. Aplicaciones de inspección y microscopía 


Robótica 


Tiempo real 

Aplicaciones en medicina, mezcla de imágenes 
inspección 

Tiempo real 

Aplicaciones en medicina, mezcla de imágenes, 
inspección 

Tiempo real 

Amplía el GVDO2 y el GM02 

Aplic. en medicina, mezcla de imágenes, inspección 


Tiempo real 
Programas en FORTRAN y BASIC 
Aplicaciones en medic., topograf., análisis 


7 seg. 


Matriz lineal 


17 seg. Inspección y microscopía, medicina 


: JU 
E S 


Fig. 7b. Continuación del cuadro comparativo de la figura 6a. 
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rio puede llevar a cabo tareas complejas a un coste moderado, sin preci- 
sar el concurso de un computador central. Sin embargo, las estaciones de- 
dicadas al procesamiento de imágenes no suelen aparecer autónomas 
(stand-alone), sino conectadas por medio de una red local como la Ether- 
net. De este modo, los usuarios pueden compartir recursos comunes a to- 
dos los usuarios de la red, como son los dispositivos de almacenamiento 
de alta capacidad (file servers), así como dispositivos de salida comunes. 
Ejemplos de este tipo de estaciones lo constituyen la estación SUN y APO- 
LLO. 

e Computadores personales. Los sistemas de procesamiento de imáge- 
nes basados en computadores personales son adecuados en aquellas apli- 
caciones en donde el bajo coste de la instalación sea un punto fundamen- 
tal. Presentan la gran ventaja de disponer de un número muy amplio de 
tarjetas de adquisición de imágenes y de software asociado. Como ejem- 
plos de este tipo de tarjetas cabe citar las de DATA TRANSLATION y las 
de PC-VISION. 


Características de los algoritmos para procesamiento de 
imágenes 


JJ 


Los algoritmos para procesamiento de imágenes deben presentar una 
característica común a todos ellos, como es su naturaleza no excesivamen- 
te compleja, puesto que deben ser aplicados a millones de datos. Como an- 
teriormente se indicó, los algoritmos para procesamiento de imágenes pue- 
den ser divididos en dos grandes grupos (figura 8): 


e Algoritmos de procesamiento puntual. En este tipo de algoritmos, las 
operaciones se llevan a cabo sobre cada pixel individualizado, sin consi- 
derar el valor de los pixels próximos. Como ejemplos cabe citar las ope- 
raciones aritméticas y booleanas entre imágenes, las transformaciones ra- 
diométricas y las operaciones estadísticas. Las transformaciones radiomé- 
tricas modifican los niveles de gris de una imagen por medio de la aplica- 
ción de funciones lineales y no lineales (incluyendo funciones aleatorias), 
dentro del rango de precisión del hardware disponible. 

Constituyen la base de una gran parte de las operaciones de procesa- 
miento, incluyendo la modificación del contraste. Las operaciones radio- 
métricas se llevan a cabo utilizando las denominadas tablas de consulta 
(look-up tables: LUT), que se encuentran situadas entre la memoria de re- 
fresco y el monitor. Las tablas de consulta se encuentran implementadas 
en memoria RAM de acceso rápido (menor de 100 nseg.), empleándose la 
intensidad de los pixels de la memoria de refresco como una dirección en 
la tabla, proporcionando un valor de salida que se presenta en el monitor. 


DD 


TRANSFORMACION 


Imagen original Imagen procesada 


(a) Transformación puntual 


TRANSFORMACION 


Imagen original “Imagen procesada 


(b) Transformación regional 


Fig. 8. Los dos tipos generales de procesamiento de imágenes. 


Este proceso se realiza a frecuencia de vídeo, permitiendo que los valores 
almacenados en la memoria de refresco puedan ser transformados en otros 
valores y representados en la pantalla de acuerdo con la función que se en- 
cuentre cargada en la tabla. Así, por ejemplo, la escala de gris puede ser 
invertida, o efectuar una transformación radiométrica logarítmica con ob- 
jeto de aumentar el contraste de las zonas. 

e Algoritmos de procesamiento espacial. En este tipo de algoritmos, la 
imagen resultante de las operaciones llevadas a cabo es función tanto de 
las intensidades de los pixels de entrada como de la de sus vecinos. 

Los algoritmos de procesamiento puntual son de fácil implementación 
en los computadores convencionales, mediante la utilización de operacio- 
nes de entrada/salida. Sin embargo, los algoritmos de procesamiento es- 
pacial requieren un mayor grado de complejidad que los anteriormente se- 
ñalados. Como ejemplos de operaciones de procesamiento espacial se en- 
cuentran las operaciones de convolución, en donde el resultado se obtie- 
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(a) Proyección de una ventana sobre una imagen que está siendo procesada. 


a pl 


Pixel a pixel dentro de una línea Línea a línea 


(b) Cálculo de ventanas consecutivas. 


Fig. 9. Ejemplos de operaciones de convolución. 


ne multiplicando cada pixel de la imagen original por un núcleo o venta- 
na de N x N pixels (figura 9). En la práctica, se utilizan núcleos de 3 x 3 
ode 5 x 5 en las convoluciones en el dominio espacial, con objeto de que 
el tiempo necesario para llevar a cabo este tipo de operaciones no sea ex- 
cesivamente grande. , 

Por otro lado, y por razones de coste, las Operaciones de procesamien- 
to espacial se realizan de forma iterativa, procesando un punto del núcleo 
o ventana por trama, con lo que se consigue un tiempo de ejecución próxi- 
mo al tiempo real. 


DL 


Un aspecto importante a tener en cuenta a la hora de evaluar los al- 
goritmos para procesamiento de imágenes es el de su tiempo de ejecución. 
Podemos considerar tres tipos de intervalos temporales: 


a) Tiempo real. Se considera que un algoritmo para procesamiento de 
imágenes se ejecuta en tiempo real si dicho tiempo es comparable al ne- 
cesario para efectuar el refresco de las memorias de imagen, de tal mane- 
ra que el usuario apenas aprecie que ha transcurrido un tiempo desde que 
lanzó la orden de ejecución y el momento en que se presenta el resultado. 

b) Próximo al tiempo real. Se considera que un algoritmo para proce- 
samiento de imágenes opera en un tiempo próximo al tiempo real si el 
tiempo de ejecución es del orden de 3 segundos. Como ejemplo de este 
tipo de algoritmos cabe señalar las convoluciones espaciales. 

c) No tiempo real. Se considera que un algoritmo para procesamiento 
de imágenes opera en un tiempo muy superior al tiempo real si el tiempo 
de ejecución está por encima de los 3 segundos anteriormente señalados. 
Este tipo de algoritmos que suponen un tiempo elevado de ejecución se 
suelen introducir en una cola de procesos. 


SISTEMAS DE PROCESAMIENTO DE IMAGENES 
BASADOS EN UN COMPUTADOR PERSONAL 


Durante los últimos cinco años, los sistemas de procesamiento de imáge- 
nes basados en ordenadores personales han comenzado a jugar un papel im- 
portante en algunas áreas de aplicación, habiendo experimentado durante 
los últimos dos años un rápido crecimiento la tecnología de este tipo de sis- 
temas. El progresivo descenso de los costes tanto del hardware como del soft- 
ware convierten al ordenador personal en una herramienta muy adecuada 
cuando el coste reducido del sistema es un requerimiento básico. 

Una configuración típica de un sistema de procesamiento de imágenes 
basado en un ordenador personal IBM PC-AT con 512 Kbytes de memoria 
principal, 20 Mbytes de memoria en disco, requeriría el siguiente hardwa- 
re adicional: 


UI 


— Tarjeta de adquisición de imágenes. 
— Cámara de vídeo y unidad de visualización. 
— Un ratón conectado vía serie. 


Opcionalmente puede añadirse el siguiente hardware: 


— Tarjeta coprocesadora conectada a la tarjeta de adquisición de imá- 
genes. 

— Tarjeta adaptadora a un bus externo, multibús, AT, etc. 

— Tarjeta de conexión a una red local: Ethernet, PC-NET, etc. 

— Tarjeta de interface a un disco óptico. etc. 
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IBM PC/AT BUS 


Fig. 10. Estructura de un sistema de procesamiento de imágenes basado en un or- 
denador personal. 
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En la figura 10 se presenta la estructura de un sistema de tratamiento 
de imágenes basado en un ordenador personal. En cuanto al software ne- 
cesario, cabe señalar que el mayor inconveniente que presenta este tipo 
de configuraciones es la limitada capacidad de direccionamiento de me- 
moria de la versión del sistema operativo MS-DOS 3.2. Este problema pue- 
de solucionarse utilizando el sistema operativo XENIX o bien en futuras 
versiones del MS-DOS. 

Una de las mayores ventajas de la utilización de este tipo de sistemas 
reside en el hecho de la gran cantidad de hardware y de software disponi- 
bles que facilitan la implementación de un sistema de tratamiento de imá- 
genes para un coste reducido. 
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APENDICE 4 


ALGORITMO DE LA TRANSFORMADA RAPIDA 
DE FOURIER 


IN 


Consideremos en primer lugar la expresión de la TDF: 


N-1 N-1 
F(u) = 1/N Y, £(x) exp[-j2ux/N] = 1/N Y, £(x) Wy (A.1) 
x=0 x=0 


siendo W, = exp[-321/N] , N= 2n y N = 2M. Seguidamente, introducire- 
mos un nuevo índice K en el sumatorio, relacionado con x mediante las 
siguientes relaciones: 


x = 2K para x par 
x= 2K + 1 para x impar 


con lo que la expresión (A.1) puede escribirse como: 
F(u) = 1/2M D) £(k) W,* = 1/211/M D, £(2K) W¿0 + 
1/M Y fQK + 1) Wi) =1/2(F,.(u) + Ej, (u) Wi) (4.2) 


par impar 


En definitiva, se ha reducido el cálculo de una transformada de N pun- 
tos a dos transformadas de (N/2) puntos. Así, por ejemplo, en el caso de 
una función con 8 puntos (£(0), f(1), ..., £(7)!, habría que considerar por un 
lado las muestras pares ¡£(0), £(2), £(4), £(6)) y por otro las impares (f£(1), £(3), 
£(5), £(7)). Cada una de las transformadas de 4 puntos se calcula como dos 
transformadas de 2 puntos. Así, para calcular la FFT de (£(0), f(2), £(4), f(6)! 
hay que dividir este conjunto en dos: la parte par (£(0), £(4)! y la parte im- 
par ¡f(2), £(6)). Análogamente se haría en el caso del segundo conjunto de 
las cuatro muestras impares: (£(1), £(3), £(5), £(7)!, que se dividiría en (£(1), 
f(S)) y 1£(3), £(7)). Así, pues, el algoritmo operaría según el esquema de la 
figura 1, de arriba a abajo: 

Seguidamente se presenta una subrutina FORTRAN, propuesta por 
Cooley y Tuckey en 1969 que permite calcular la FFT mediante el método 
descrito anteriormente. Los parámetros de la subrutina son dos. F es un 
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£(0) 


f(1) £(2) 


f(6) £(1) £(5) £(3) £(7) 


DOS PUNTOS 


TRANSFORMADAS DE 


CUATRO PUNTOS 


TRANSFORMADAS DE 
OCHO PUNTOS 


FFT 


Fig. 1. Ordenación de los datos de la función Íf(0), f(1), ..., f(7)! para el cálculo de 


la FFT. 


array cuya transformada se desea calcular y LN = n, siendo N = 2n la di- 
mensión del array considerado. Como resultado de los cálculos esta sub- 
rutina devuelve en F el valor de la transformada de Fourier. Obsérvese 
que F ha sido dimensionado en la subrutina como un array de números 
complejos, luego si partimos de un array F correspondiente a una función 
real, la parte imaginaria de F debe ser cero antes de llamar a esta subrutina. 
SUBROUTINE FFT(F,LN) 

COMPLEX F(1024),U,W,T,CMPLX 
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N — 


PI=3.141593 
N=2**LN 
NV2=N/2 
NM1=N-1 

J=1 

DO 3 I=1,NM1 
IF(I.GE.J) GO TO 1 


IF(K.GE.J) GO TO 3 
J=J-K 


ad + 


K=K/2 

GO TO 2 

J=J+K 

DO 5 L=1,LN 
LE=2**L 

LE1=LE/2 
U=(1.0,0.0) 
W=CMPLX(COS(PI/LE1),-SIN(PI/LE1)) 
DO 5 J=1,LE1 

DO 4 I=J,N,LE 
IP=I+LE1 
T=F(1P)*U 
F(IP)=F(D)+T 
F(D=F()+T 
U=U*W 

DO 6 I=1,N 
F(D)=F()/FLOAT(N) 
RETURN 

END 
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APENDICE Y, 


PROGRAMAS DE EXPERIMENTACION 


El objetivo de los programas que se presentan en este apéndice consis- 
te en proporcionar una herramienta sencilla, flexible y fácil de utilizar por 
aquellos usuarios de ordenadores personales. Los programas se han escri- 
to en Fortran 77, requiriendo para su utilización una impresora de carac- 
teres con objeto de poder imprimir los resultados. 


Las imágenes utilizadas tienen un tamaño de 64 x 64 pixels, almace- 
nando cada pixel como un valor complejo en punto flotante. Este tamaño 
de imágenes es muy apropiado para su presentación tanto en un monitor 
alfanumérico como para su impresión en una impresora de líneas, propor- 
cionando suficiente resolución para el fin didáctico a que se dirige. Este 
formato de datos complejos se ha efectuado con objeto de facilitar el aná- 
lisis de Fourier. La precisión de punto flotante asegura la no introducción 
de errores de redondeo. Por otro lado, el cálculo de la transformada de 
Fourier se ha efectuado de tal manera que la frecuencia cero corresponde 
al centro de la imagen presentada. La presentación en la pantalla se ha he- 
cho utilizando una sola línea de caracteres para producir ocho niveles de 
seudogris. 


En el programa se han definido tres buffers de imágenes de 64 x 64 
pixels, lo cual permite la realización de operaciones de suma, escalado, 
transformación, etc. La utilización de este programa de entrenamiento es 
sencilla, puesto que se basa en la elección entre una lista de comandos pre- 
definidos acompañados de uno o dos argumentos correspondientes a los 
«buffers» de imagen. El resto de los parámetros se introducen en respues- 
ta a las peticiones efectuadas por el programa, asignándose valores por de- 
fecto en el caso de no especificar dichos parámetros. 


Asimismo, el programa permite el diagnóstico de errores de sintaxis en 
los comandos, devolviendo control al usuario para su nueva introducción. 
Los errores producidos como consecuencia de la introducción de un pa- 
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rámetro erróneo implican la necesidad de introducir de nuevo el coman- 
do. En la figura 1 se presenta una imagen de barras y en la 2 el módulo de 
su transformada de Fourier. En la 3 se presenta una imagen constituida 
por un cuadrado de niveles de gris y en la 4 el módulo de su transformada 
de Fourier. 


PROGRAM AUTOMATA 


Este programa es el que maneja por media 

de una autómata las —cperaciones que se 
pueden efectuar en el procesamiento de imagenes. 
El autómata maneja una tabla multiple que en 
ayuda de un scannner se encarga de decodificar 
una linea de comandos entregada a él. 


000000 


COMMON /BUFFER/ A,B,C 
COMPLEX—AC4096),B(40767,€ (4076) 
COMPLEX AA(64,64) ,EB(54,64),CC (64,66) 
EQUIVALENMCE CA,AAI,1B,BB>,1C0,C0C) 


CHARACTERXI COMANDO(13,30)7,LIMEACBO),MANDATOC13> 
ECHARACTER+*80-—— LIN 

CHARACTERx220 FAL 

CHARACTER*22-—PALABRASCIO) 

CHARACTER+*13  NOMBRES(30), TOKEN, NOMBRE 
EQUIVALENCE (MOMBRES COMANDO), CTOKEN, MANDATO? 
EQUIVALENCE (LIN,LINEA), (PALABRAS, PAL) 

COMMON /COMMAND/LIN, MANDATO 
COMMON-—7/COMMANDIA/LONGITUD: INDICE 

INTEGER DIRECCION(30) 

direccion que se debe ejecutar despues del comanda 
INFTEGER—REGRESOCIO? 

e direccion de retarna uña vez que se ejecuta 
INTEGER SUMA,RESTA, MULTIPLICACIÓN, DIVISION, DIRECTA, INVERSA 
INTEGER ESCRIBIR, LEER,OPERACION, SALTR,CARACTER 
INTEGER REACOIMAG,COMJU, RECOPOC,POCOREC,MEG 
INTEGER LOGAE,LOGALO,REALO, IMAGO 
INTEGER CERO, CAMBIO, CUADRADO, ORGO,DFFT, INVFFT 


E Datas para las constantes que se ocuparan en el programa. 


DATA -SUMA,RESTA,MUETIPLICACION, DIVISION/1,2,3,47 

DATA TESCRIBIR,LEER71527 

DATA SALIR/1/ 

DATA REALOIMAG,CONJU,RECOPOL,POLOREC, 

FNEG,LOGAE ,LOGA1O,REALO, IMAGO/5,6,7,8,9,10,11,12,13/ 

DATA CERO, CAMBIO, CUADRADO, ORGO,DFFT,INVFFT/14,15,16,17,18,197 


= Máximo númerc de cemandos que tiene la tabla 
DATA MAXIND /10/ 
£ PMI III 


£ Comandas que son aceptados. 
DATA NOMBRES/*DISCO y 
* *SUMO => 
* “RESTO == 
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"RESTAR = 
"DIVIDIR E 
"DIVIDO = 
> MULTIPLICO - 
MULTIPLICAR >, 
"NEGAR = 
"SALIRSE = 
*SUME = 
"SUMAR = 
—RESTE == 
"DIVIDE 5 
"MULTIPLIQUE>, 
"NOT = 
"NIEGUE a 
"TERMINAR = 
> ACABAR 7 
"EXIT = 
"END = 
"LOGOUT z 
"CAMBIE = 
*TRASLADE ; 
"TRASLADO : 

. 

, 

hd 

, 

7 


CARRETAS 


"CAMBIAR - 
"LOGARITMO : 
"CERO: - 
“CLEAR : 
"SORT z 


FA A RA AAA A A A RA A AR AAA AAA 


E Palabras que se utilizan en una linea de comandos. 
DATA PALABRAS 

* "BASE -10 REAL-IMAG-R-1-” 
*”REAL— IMAGINARIO COMJUG” 
*”ADO-REC=POL RECTANGULA”? 
*"R POLAR RP PR 2 ORIGEN” 
** ORO DER=IZ0 IZQUIERDA” 
** DERECHA TRASLADE 
> NEGAR PARTE BU” 
EFFER— EERESCRIBIR 
$ 
e 


» 


. 


Quo ss ss 


l Direcciones del regresa. 

DATA REGRESO/24,24,24,24,24,24,24, 
*24,24,24,24,26,24,24,24,24,24,24, 
*24,24,24,24,24,24,24,24,24,24,24,247 

€ Direcciones dende voy despues de decodificar el autamata. 

DATA DIRECCION 
A O O 
*3,6,47,3,13,13,1,1,1,1, 

A A A A 


€ COMIENZO” EL PROCESO CON EL AUTOMATA. 
30 CONTINUE . 
CALL CLRSCR 


2) HAHAHAHA PRUEBA ERARIO 
DUTISR=EG3 


ral 


20 


10 


DOOoOn+ 


CALE-PATRON(KA) 
DO 13 I=0,11 
DO-14=J=1,64 
PRINT *,REAL(A( I*b4+J)) - 
CONTINUE 
PRINT-+* 
CONTINUE 
CONTIMUE 
PRINF- FIN? 
sTorP 
IN II IA 
CALL LEERLINE 
INDICE=1 
CALL SCAMMER 
JUMF=1 
J=1 
CONTINUE 
LUGAR=1 
k=ICHAR(COMANDO(1,3)) 
L=ICHARIMANDATO(LUGARI) 
TE AG ES) FHEN 
Veremos si este es el comando que se invoca 
1=2 
CONTINUE 
E=ICHARIMANDATO(TI)) 
L=ICHAR(COMANDO(1,J)>>) 
IF (K.EQ.L) LUGAR=LUGAR+1 
I=I+ 
IF (LUGAR.EQ.LONGITUD) JUMP=DIRECCION(J) 
TF ALUGAR-EO.CONGITUD)—GO-TO-—40 
IF (I.LE.LONGITUD+1) GO TO 10 
ENDIF 
J=J+1 
IF (J.LE.MAXIND) GO TO 20 
G0-TO-100 
CONTINUE 


A A A A A A A A A 
Comienzo el procesamiento del comando dada. 


1F ((JUMP.LT.1).AND. (JUMP.GT.24))60 TO 100 
G0-TO-(100,200,300,400,500,600,700,800,900,1000,1100, 


*1200,1300,1400,1500,1600,1700,1800,1900,2000,2100, 
*2200,2300,30-), JUMP 
0 0 0 o JD JE A e JJ EI JJ MJ Jl Jl JE JE O E JE JE E JE E E E e 


E Aqui nunca podre regresar de la salida del automata.  * 
pe * 
ec * 


PIMIENTO III III III 


100 CONTINUE 


OPERACION=5 
CALL FILES(OPERACION, NOMBRE, A) 
PRINT *,LIN 
PRINT-x 
G0 TO 30 
c FMI MNR DIM III DD DD MD TI EH A A ll le Bl Dl RA 


200 CONTINUE 
300 CONTINUE 


TZ 


o0 


CONTINUE 

CONTINUE 

Los comandos invocados son 206,300,400 y 500 para suma, 

resta,multiplicación y división. 

CALL SCANNER 

CARACTER=ICHARIMANDATOC1)) 

PRINT *,MANDATO(1),TOFEN, INDICE 

PAUSE "ESTOY EN 200 AL S00” 

IF ACARACTER.EO.65)—THEN 

Era el buffer A. 
CALL SCANNER 
PRINT-*,MANDATO (1), INDICE, TOFEN,1ichartmandateat1)) 
PAUSE ”* INGRESO AL TIPO DE OPERACION” 
OPERACION=SALIR 
CARACTER=ICHAR(MANDATO(1)) 
IF (CARACTER.EO.43) DPERACION=SUMA 
TF <CARACTERSEO 45) OPERACION=RESTA 
IF (CARACTER.EQ.42) OPERACTON=MULTIFLICACION 
TF 1CARACTER.£0.538) OPERACION=DIVISION 
IF (CARACTER.EO.47) OFERACION=DIVISION 
IF COPERACION.E0Q.SALIR) GO TO-100 

Erg un signo de operacion válido: 

PRINT-X,INIDICE 
CALL SCANNER 
CARACTER=ICHAR(MANDATO(1)) 
PRINT-*,MANDATOC1), INDICE, TOFEN, OPERACION 
PAUSE *DECODIFTIOUE-LA- OPERACION VEREMOS DOMDE VAMOS” 
IF (CARACTER.EO.65) CALL FUNCION(OFPERACION,A,A) 
IF (CARACTER.E0.66) CALLE FUNCIONCOPERACION, AE) 
IF (CARACTER.EQ0.67) CALL FUNCION(OPERACTON,A,C) 
IF ICCARACTER.NE-.65)7.0K.(CARACTER.NE.66).0R.: 


hd (ECARACTER.NE.67)) 60 TO 1060 


JUMP=REGRESOCIUMP> 
PRINT *, JUMP 
PAUSE” "VEREMOS DONDE-SACTO” 
GO TO 46 

ENDIF 

IF (CARACTER.EQ0.66) THEN 

Era el buffer B. 
CALL SCANNER 
OPERACION=SALIR 
CARACTER=ICHAR(MANDATO (13) 
IF (CARACTER.E0.43) OPERACION=SUMA 
IF (CARACTER.EQ.45) OPERACION=RESTA 
IF (CARACTER.E0.42) OPERACION=MULTIPLICACION 
IF (CARACTER.EO.58) OPERACION=DIVISION 
IF (CARACTER.E0.47) DOPERACION=DIVISION 
1F (OPERACION.EGO.SALIR) 60 TO 100 


Eracun signo de cperacion válido. 
CALL SCANNER 
CARACTER=ICHARIMANDATO(1)) 
TF (CARACTER: EO.65)- CALL FUNCIONCOPERACION, EA) 
IF (CARACTER.EO.66) CALL FUNCION(OPERACION,B,E) 
TF (CARACTER.EO.67) CALC FUNCIONCOPERACION, EC) 
IF (¿CARACTER .NE.65).OR. (CARACTER.NE.66).0R. 

(CARACTER NE=26777—G0-TO-100 


JUMP=REGRESO (JUMP) 
60-TO—40 


TS 


ENDIF 
IF ACARACTEREO677—THEN 
6 Era el buffer E. 
CALISCANMNER z 
OPERACTION=SALTR 
CARACTER=ICHAR(MANDATO(1)) 
TF ACARACTER.E0.43) OPERACION=SUMA 
IF (CARACTER.EO.45) OPERACION=RESTA 
IF (CARACTER.E£Q0.42) OPERACTON=MULTIFLICACION 
IF ACARACTERSEC5 59) OPERACTON=DIVISION 
IF (CARACTER.EQ.47) OPERACION=DIVISION 
TF COPERACIONS ES: SACITR)—=G0-TO-100 
E Era un signo de aperacian válido. 
CALL SCANNER 
CARACTER=ICHAR<MANDATO 41) 
IF (CARACTER.EQ.6353) CALL FUNCION(OPERACION,C,A) 
TF (CARACTER; EQ, 66) CALL FUNCIONO PERACION,C,B> 
IF (CARACTER.EQ.67) CALL FUNCTON(OPERACTON,C,C) 
IF IACARACTERME-.63).0R.1CARACTER.NE.66)7.0K. 


* (CARACTER-.NE.67)) 60 TO 100 
JUMP=REGRESOCIUMP> 
GO TO 40 
ENDIF 
GO TO 100 
-— MIN N III III MITIN INICIO 
e 
600 CONTINUE 
£ Comando para teercen—buffers, 
CALL SCANNER 
OPERACION=LEER 
1F-(LONGIFUD.6F.1)-60-TO-100 
IF (CARACTER.EQ.65) THEN 
NOMBRE="EUFFERASDAT” 
CALL FILES(OFERACION, NOMBRE,A) 
JIUMP=REGRESOCIUMP) 
60 TO 40 
ENDIF 
IF (CARACTER,EQ.66) THEN 
NOMERE=" BUFFERE.DAT” 
CALL FILES(OPERACION, NOMBRE,E) 
JUMP=REGRESOCIUMP) 
GO TO 40 
ENDIF 
IF (CARACTER+.E0,67) THEN 
NOMBRE=”"BUFFERC.DAT” 
CALLE FTICESCOPERACION, NOMBRE, C) 
JUMP=REGRESO (JUMP) 
G0-TO-40 
ENDIF 
ce Errar en lectura de buffers. 
GO-TO-100 
e] 
= Comando para escribir las buffers en disca. 


700 CONTINUE 
CALL SCANNER 
OPERACION=ESCRIBIR 
IF (LONGITUD.GT.1) GO TO 100 
IF <CARACTER.EO0 65) THEN 
NOMBRE=" RUFFERA.DAT”* 
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CALL FILES(OPERACION, NOMBRE, A) 
JUMP=REGRESO( JUMP) 
60 TO 40 
ENDIF 

IF (CARACTER.EQ.66) THEN 
NOMERE="EUFFERB.DAT* 
CALL —FILESCOPERACION, NOMBRE, El) 
JUMP=REGRESO (JUMP) 
GOTO 40 
ENDIF 

IF (CARACTER.EQ.67) THEN 
NOMBRE="BUFFERC.DAT" 
CALL FILES(OPERACION,NOMBRE,C) 
JUMP=REGRESOCIUMP) 
60 TO 40 
ENDIF 


A — continuación estan todas las funcicnes que es posible 
hacer con laos buffers, SUMA, RESTA, DIVISION y MULTIPLICACION. 
CONTINUE 
CALL SCANNER 
IF (CARACTER.E0.65) THEN 
CALL SCANNER 
OFERACTON=SAL.IR 
IF (CARACTER.EQ.43) OPERACION=SUMA 
IF (CARACTER.EQ.45) OPERACITON=RESTA 
IF (CARACTER.E0.42) OPERACION=MULTIPLICACIOÓN 
IF (CARACTER.EO.58) OPERACION=DIVISION 
IF (CARACTER.EQ.47) OPERACION=DIVISION 
IF (OPERACION.EO.SALIR)-—60- TO 100 
CALL SCANNER 
IF (CARACTER.EO.65) CALL FUNCIONCOPERACION,A,A) 
IF (CARACTER.EQ.66) CALL FUNCION(OPERACION,A,H) 
IF (CARACTER.E0.67) CALL FUNCIONCOPERACION,A,C) 
IF ((CARACTER.ME.65) OR. (CARACTER.NE.66) .OR. 
(CARACTER.NE.67))—G0-—TO-100 
JUMP=REGRESO( JUMP) 
GOTO 40 
ENDIF 
IF (CARACTER.EOQ.665) THEN 
CALE-SCANNER 
OFERACTON=SAL TR 
IF (CARACTER.E0.43) OPERACION=SUMA 
IF (CARACTER.EO.45) OPERACION=RESTA 
IF (CARACTER.EQ.42) OPERACTON=MULTIFLICACION 
TF (CARACTER, E0O.58) OPERACION=DIVISION 
IF (CARACTER.EO.47) OPERACION=DIVISION 
IF COPERACION.EO. SALIR) 60 TO-100 
CALL SCANNER 
IF (CARACTER.EO.65) CALL FUNCIONTOPERACION, E, A) 
TF A CARACTER, E0O.66) CALLE FUNCIONCOPERACION, E, E) 
IF (CARACTER.EGQ.67) CALL FUNCION(OPERACION, R,C) 
TF ACCARACTER MES 6370. (CARACTER.ME-66).0K. 
(CARACTER.NE.67)) 60 TO 100 
JUMP=REGRESOCIUMP) 
60 TO 40 
ENDIF 
TF A CARACTER EO.67)—THEN 
CALL SCANNER 
OPERACTION=SALTR 
IF (CARACTER.EO.43) OPERACION=SUMA 


TS 
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IF (CARAC 


TER.EO0.45) OPERACTON=RESTA 


IF (CARACTER. E0O. 42) OPERACION=MULTIPLICACION 


IF (CARAC 


TER.EGO.58) OPERACION=DIVISION 


IF (CARACTER.E0.47) OPERACION=DIVISION 


IF-COPERA 


CIONEOSALTR)—60-TO-100 


CALL SCANNER 


IF (CARAC 
IF-ACARAC 
IF (CARAC 


TER.E0.65) CALL FUNCIONCOPERACION,C,A) 
FER +E0.66) CALL FUNCIONCOPERACION, EE) 
TER.EO0.67) CALL FUNCION(OFERACION,C,C) 


TF AACARACTER<NE-6570RCEARACTERNE-667.085 
(CARACTER.NE.67)) GO TO 100 


JUMP=REGRI 
60-TO-40 
ENDIF 


ESOCIUMP) 


Procedimiento para obtener filtro en imagen c bien FFT directa 


inversa. Ademas cambics varics en una matriz. 


Intercambio parte REAL par IMAGINARIA 


CONTINUE 
OPERACION=REALOIMAG 
CALL SCANNER 
GO0-TO-1 


Procedimiento para cbtener el conjugado de una matriz. 


CONTINUE 
OPERACION=CONJU 
GO TO 1 


Procedimiento para 
CONTINUE 
DPERACION=RECOPOL 
60 TO 1 


Procedimento para 
CONTINUE 
OPERACION=POLOREC 
GO TO 1 


Procedimiento para 
CONTINUE 
OPERACION=NEG 

Go TO 1 


Procedimento para 

Procedimiento para 

iz 

CONTINUE 

CONTINUE 

CALL SCANNER 

K=INDEX(FAL, TOKEN) 

TEIEZEOIITHEN 
OPERACI 


cbtener el cambio de rectangular a polar: 


convertir de polar a rectangular. 


negar una matriz. 


sacar el togarítmo en base £ de una matriz 
cbtener el logaritmo en base 10 de una matr 


ON=SALIR 


CALL SCANNER 


k=INDEX 


PAL TOFREN> 


IF (K.E0.4) OFERACTON=LOGAE 


1900 


IF (K.E0.6) OPERACION=LOGA1O 
IF(OPERACION.EQ.SALIR) GO TO 100 
GU-—FO-—1 
ELSEIF(K.E0.4) THEN 
OPERACTON=LOGAE 
GO TO 1 
ECSETEUEOTS7THEN 
OPERACTON=LOGA1O 
GO TO 1 
ENDIF 
GO TO 100 


Procedimiento para llevar a cerco la parte REAL de una matriz. 

o la parte IMAG 

CONTINUE 

CONTINUE 

CALL SCANNER 

K=INDEX(PAL,TOKEN) 

IF(K.EO.150) THEN 
CALL SCANNER 
OPERACION=SAL IR 
K=INDEX(PAC, TOKEN) 
IF(K.EO.11) OPERACION=REALO 
IF(K.E0. 14) OPERACION=IMAGO 
IF(SALIR.EQ.OFERACION)GO TO 100 
GOTO 

ELSETFA EO-11)—THEN 
OPERACION=REALO 
G0-TO-1 

ELSEIF(K.E0.14) THEN 
OPERACION=IMAGO 
60 TO 1 

ENDIF 


Procedimiento para borrar el buffer indicado. 

CONTIMUE 

CALL SCANNER 

K=INDEX(PAL, TOKEN) 

IF(kE0.58)—THEN 
OPERACION=CERO 
GO-TO1 

ELSEIF(K.E0.153) THEN 
GO TO 1 

ELSEITF(K.E0O.1) THEN 
INDICE=INDICE-1 
GO-TO-1 

ELSETFKk.E0Q.2)-—THEN 
INDICE=INDICE=1 
60-—TO-—1 

ELSEIF(K.EQ.39) THEN 
INDICE=INDICE—1 
60 TO 1 

ENDIF 

G0-TO-100 


Procedimiento para intercambiar parte derecha por 
la izquierda en una imagen. 

CONTINUE 

OPERACTON=CAMBIO 

GO TO 1 


AL 


ES 


0000000000 


Precedimiente para sacar el cuadrado de la magnitud. 
CONTINUE 

OPERACION=CUADRADO 

GO TO 1 


Procedimiento para poner en cera el origen de la imagen. 
CONTINUE 

OPERACION=0RGO 

GO0-TO-1 


Pracedimienta para cbtener la transformada 
de FOURIER de la imagen 

CONTINUE 

OPERACION=DFFT 

GO TO i 


Procedimiento para obtener la tranformada inversa de FOURIER. 


CONTINUE 
OPERACITON=INVFFT 
Go TO 1 


Pracedimiento para procesar todos las comandos indicados. 


CONTINUE 

CALL SCANNER 

IF (LONGITUD.GT.1) 60 TO 100 

IF (CARACTER.EQ.£53) CALL FUNCION(OPERACION,A,B) 
IF <CARACTER:EO.66) CALL FUNCIONO PERACTON, EA? 
IF (CARACTER.EQ.67) CALL FUNCIONCOPERACION,C,A) 
IF IACARACTER.NE-657.0K.(CARACTER.NE.66)7.0K: 
* (CARACTER.NE.67)) THEN 


GO-TO-100 
ELSE 
JUMP=REGRESO (JUMP) 
GOTO 40 
ENDIF 
JUMP=1 
G0-TO-—40 
END 
END 


SUBROUTINE FUNCION(CODIGO,A,B) 
ESTA SUBROUTINE EFECTUA LAS OFERACINES ARITMETICAS 
SOBRE EL PLANO IMAGEN DE MATRICES A Y B AMBAS COMPLEJAS. 


CODIGO=1 
SE-EFECTUA-TA-—SUMA-DE- MATRICES A=A-+-B 


copIco=2 
SE -EFECTUA-LA-RESTA DE MAPFRICES A=A-——E 


2000000000000000000000000000000000000000000000000O 


10 


CcOoDIGO=3 
SE EFECTUA LA MULTIPLICACION DE MATRICES A=A * E 


CODIGO=4 
SE EFECTUA LA DIVISION DE MATRICES A=A / B 


CODIGO=5 
SE EFECTUA EL” TRASCADO DE A ==> B YB ==>A 


CODIGO=6 
SE-EFECTUA EL” TRASCADO DE-PARTE-REAL— POR—IMAGINARIA 


CoDIGO=7 
SE EFECTUA EL CAMBIO DE - MATRIZ POR EL CONJUGADO: 


CoDIGO=8 
SE EFECTUA EL CAMBIO DE RECTANGULAR A POLAR. 


CODIGO=9 
SE EFECTUA EL CAMBIO DE FOLAR A RECTANGULAR. 


CODIGO=10 
SE-EFECTUA EL NEGADO” DE LA MATRIZ 


CODIGO=11 
SE-EFECTUA EL LOGARITMO EN BASE-E: 


CoODIGO=12 
SE EFECTUA EL LOGARITMO EN BASE 10. 


CODIGO=13 
SE EFECTUA LA PUESTA A CERO DE LA PARTE REAL. 


CODIGO=14 
SE-EFECTUA CA PUESTA A CERO DE LA PARTE IMAGINARIA 


CODIGO=13 
SE BORRA EL BUFFER. 


CODIGO=16 
SE INTERCAMBIA PARTE DERECHA POR IZQUIERDA EN LA IMAGEN. 


CODIGO=17 
SE SACA EL CUADRADO DE LA MAGNITUD. 


CODIGO=18 
SE-PONE-—A- CERO EC ORIGEN: 


COMMON —7COMMAND/CINEA, COMANDO 

COMPLEX RESP(4096),AM(64,64) ,EB(64,64) 
COMPLEX A(4096),B(4096) 

EQUIVALENCE (RESP ,AA,BB) 

INTEGER CODIGO 


IF(CODIGO.EQ.1) FHEN 

SUMA-DE- MATRICES COMPLEJAS 
DO-10—1=1 409651 
ALI)=A(1I) + BI) 
CONTIMUE 


TO 


ELSETF”(CODIGO.E£EC0.2)— THEN 


E RESTA DE MATRICES 
DO-20-1=1,40965,1 
ACI)=ACI) — BClI) 5 
20 CONTINUE 
ELSEIF (CODIGO.E0.3) THEN 
E MULTIPLICACION DE MATRICES 
DO 30 1=1,4096,1 
ACI)=ACI) * BC1) 
30 CONTINUE 
ELSEIF (CODIGO.EQ.4) THEN 
E DIVISION DE MATRICES 
DO-40—I=1,4096,1 
ALI)=ACI) / BL) 
40 CONTINUE 
ELSEIF (CODIGO.EQ.5) THEN 
a SE-HACE EL INTERCAMBIO DE A --> BY B-=->4 
DO-50-—1=1,4096,1 
RESP(1)=A11) 
ACD=BCD) 
B(1)=RESP(1) 
50 CONTINUE 
RETURN 
ELSEIF (CODIGO.EO.6) THEN 
E INTERCAMBIO PARTE-REAL POR IMAGINARIA 


PRINT-*,"—Intercambic de parte REAL por IMAGINARIA ” 
DO 101 1=1,4096 
PREAL=REALCACIO) 
FIMAG=AIMAG(A(I)) 
ACI)=CMPLX(PIMAG,PREAL) 
101 CONTINUE 


RETURN 
ELSEIF (CODIGO.E0.7)— THEN 
e CAMBIO LA MATRIZ POR EL CONJUGADO. 


200 PRINTx," cambic matriz por el CONJUGADO” 
DO 201 1=1,4096 
ACI)=CONJG(A(I)) 

201 CONTINUE 


RETURN 
ELSEIF (CODIGO.E0.8) THEN 
e CAMBIO DE RECTANGULAR—A-POLAR:. 


300 PRINT x,” Cambio de rectangular a palar ” 
DO-301—1=1,4096 
VALOR=CABS(A(1)) 
ANGULO=90.0XSIGN(1.O,AIMAGCACIT)) 
PREAL=REALCACI)) 
PIMAG=AIMAG(A(I)) 
IF (RESBCCA CI) NE-0.0)—ANGUELO=57.29578*ATANSPIMAG;PREAL) 
A(1I)=CMPLX (VALOR, ANGULO) 
301 CONTINUE 


RETURN 
ELSEIF (CODIGO.EQ.9) THEN 
ce CAMBIO DE POLAR A RECTANGULAR: 


400 PRINT x*«,7 Convierto de polar a rectangular 
DO-401—1=1,4096 
ANGULO=0.01745329*AIMAG(A(I)) 
PREAL=COSTANGULO)*REALCACI)) 
PIMAG=SIN(ANGULO)*REAL(A(1)) 
A(1)=CMPLX(PREAL,PIMAG) 

401 CONTINUE 


30 


2002 


RETURN 


ELSEIF (CODIGO.EQO.10) THEn 


MEGADO DE LA MATRIZ. 


PRINT*,>El negado de la matriz” 


DO-501-—T=1,4096 
ACI)=-A(I) 
CONTINUE 


RETURN 


ELSEIF (CODIGO.EC.11)— THEM 
LOGARITMO EN BASE E 


PRINT Xx,” Lagarítmo en base E ” 


DO—601-—T=1,4096 

VALOR=REAL(A(I)) 

IF 1WACOR.E0.0.0) VACOR=10.0E=20 
VALOR=ALOG(VALOR) 
ACI)=CMPLX(VALOR,O.0) 

CONTINUE 

RETURN 


ELSETF (CODIGO.EO.12)—THEN 


LOGARITMO EN BASE 10. 


PRINT-%,* LCogaritmo en base 10” 


DO 701 T=1,4096 

VALOR=REAL(A(I)) 

IF (VACOR.E0.0.0) VACOR=10.0E-20 
VALOR=2 . 30259*ALOG(VALOR>) 
ACI)=CMPLX(VALOR,O.O) 


CONTINUE 


RETURN 


ELSEIF —(CODIGO.E£0.13)— THEN 


PONER A CERO LA PARTE REAL. 


PRINT-*," Poner a cerco la parte real” 
DO 801 I=1,40965 

VALOR=AIMAGCACI)) 

ACI)=CMPLX(0.0, VALOR) 

CONTINUE 

RETURN 

ELSEIF (CODIGO.EO.14) THEN 
PONER—A-CERO- LA PARTE-—IMAGINARIA: 
FRINT x,* Poner a cera la parte imaginaria” 
DO—901—T=1,4096 

VALOR=AIMAG(A(I)) 
ACTI=CMPEXIVACOR, 0-0) 


CONTINUE 


RETURN 
ELSETF (CODIGO. EO.15)7—THEN 
BORRAR EL BUFFER. 


PRINT*, Barro el buffer >” 


re 


DO 1001 I=1,4096 

ACI=CMPLX(0.0,0.0) 

CONTINUE 

RETURN 

ELSETF”(CODIGO.EO.16)7— THEM 
INTERCAMBIO PARTE DERACHA POR IZOUIERDA DE LA IMAGEN. 
PRINT-*X,7 Intercambio parte derecha”, 
per—parte izquierda en la imagen” 
DO 2002 I1=1,4096 

RESPCI=ACI) 

CONTINUE 

DO 2001—1=1,32 


81 


3000 


3001 


nonanidaia 


10 


82 


11=65-1 

DO 2001 J=1,64 

BAR=AACIZIO) 

AA(T,1)=AA(J,11) 

AACI,II=BAK 

CONTINUE 

DO 2003 I=1,4096 

ACI)=RESPCTD) 

CONTINUE 

RETURN 

ELSEIF- (CODIGO.E0.17)-—THEN 
CUADRADO DE LA MAGNITUD. 

PRINT-*," Cuadrado de la magnitud ” 
DO 3001 I=1,4096 
ACTI=CABSTACIA*2 

CONTINUE 

RETURN 

ELSETF-(CODIGO.E0O.18)—THEN 

PONER A CERO EL ORIGEN DE LA IMAGEN. 
PRINT>*,* Poner a cero el crigen > 
A(1089)=10.0,0.0) 

RETURN 

ENDIF 

END 

SUBROUTINE-PATROMtPAT,A) 


ESTA RUTINA GENERA EN LA MATRIZ *A” UN PATRON DE TEST: 


PAT==1 ==5 FRAMIAS HORIZONTALES 
FAT = 2 ==> FRANJAS VERTICALES 
FAT = 3 ==> AJEDREZ 


INTEGER FAT,CONTA 
REAL BLANCO 

REAL NEGRO 

REAL VALOR, IMA 
COMPLEX—A(4096) 


BLANCO =-10.0 
NEGRO = -10.0 
IF CPAT-EO-1) THEN 
PATRON-DE FRANJAS HORIZONTALES 
CONTA=0 
YVACOR=ECANCO 
DO 10 I=1,4096,1 
CONTA=CONTA+1 
1IF-(CONTA.EO.256) THEN 
CONTA=0 
VALOR = =VALOR 
ENDIF 
IMA=SAIMAGCACI)) 
A(CI)=CMPLX(VALOR, IMA) 
CONTINUE 
ELSEIF (PAT.E0.2) THEN 
PATRON DE FRANJAS VERTICALES 
CONTA=0 
VALOR=ELANCO 


DO -20-—1=1,4076,1 
CONTA=CONTA+1 
IF (CONTA.EQ.4) THEN 
CONTA=0 
VALOR = VALOR 
ENDIF 
IMA=AIMAG(A(I)) 
ACTIr=CMPEXIVACOR, IMA) 


20 CONT INUE 
ELSEIF (PAT.EQ.3) THEN 
€ PATRON DE AJEDREZ 


DO 30 I=1,4096,1 
TEZOCIEMOD 27 EOI THEN 
IF-—CCFILA MOD 2),E0.1) THEN 
VALOR = BLANCO 
ELSE 
VALOR = NEGRO 
ENDIF 
ELSEIF ((FILA MOD 2).E0.1) THEN 
VALOR = NEGRO 


ELSE 
VALOR = BLANCO 
ENDIF 
A(T)=CMPLX (VALOR, VALOR > 
30 CONTINUE 
ENDIF 
RETURN 
END 
c FEMME MMM MM MMMM MM MMM 


SUBROUTINE FILES(OPERACION, ARCHIVO, MATRIZ) 


ESTA RUTINA HACE-LA APERTURA, CERRADO, LECTURA 
O ESCRITURA DE UN ARCHIVO. LA LECTURA O ESCRITURA 
SE HACE POR MEDIO DE LA MATRIZ E. 


OPERACION=1 
ENTONCES SE ABRE UN ARCHIVO, 


OPERACION=2 
ENTONCES SE CIERRA UN ARCHIVO. 


OPERACION=3 
ENTONCES SE LEE UN ARCHIVO. 


OPERACION=4 
ENTONCES SE ESCRIBE UN ARCHIVO. 


OPERACION=S 
ENTONCES SE ESCRIBE UN MENSAJE 
DE COMANDO MAL FORMADO. 


LOS DATOS LEIDOS O HA ESCRIBIR SE DEBEN 
ENTREGAR EN LA MATRIZ B DEL 
BUFFER ARCH. 


0000000000000 0000000000000 


» 
RAEE EE ER EEK EEE RE EMIR 


COMPLEX B(4096),BB(64,64) ,MATRIZ(40965) 
EQUIVALENCE (E, BB) 


e 


e 


e 


10 


11 


21 


20 


100 


84 


CHARACTER+*12 ARCHIVO 
INTEGER OPERACION 


IF (OPERACION.EC.1) THEN 
ABRO UN ARCHIVO. 

OPEN(UNIT=4,FILE=ARCHIVO, 
*  ACCESS="DIRECT” ,RECL=512) 

RETURN 

ELSEIF (OPERACION.EOQ.2) THEN 
CERRAR UN ARCHIVO. 

CLOSE (UNIT=4) 

RETURN 

ELSEIF (OPERACION.EO.3) THEN 
SE LEE UN ARCHIVO EN LA MATRIZ DEL BUFFER ARCH LLAMADA B. 

DO—10—T1=1,64,1 

DO 10 J=1,64,1 

READ(UNIT=4,REC=I)BB(1,J) 

CONTINUE 

DO 11 I=1,4096 

MATRIZCI)=BC1) 

CONTINUE 

RETURN 

ELSEIF COPERACION.EQ.4) THEN 
ESCRIBO EL BUFFER ARCH EN EL DISCO. 

DO 21 1=1,4096 

B(1I)=MATRIZ(T) 

CONTINUE 

DO 20 I=1,64,1 

DO 20 J=1,64,1 

WRITE CUNIT=4,REC=I>EB(1,J) 

CONTINUE 

RETURN 

ELSEIF (OFERACION.EQ.5) THEN 
ESCRIBO UNA PANTALLA DICIENDO QUE EL ERROR SE 
DEBE A ILEGAL COMANDO. 

CALL CLRSCR 

PRINT-*,” ERROR: el comando”, 
*? que invoco es ilegal.” 

RETURN 

ENDIF 


RETURN 
END 


SUBROUTINE-—CERSCR 

CHARACTER*x1 FILA(10) 

CHARACTER*10 BORRAR 

EQUIVALENCE (FILA,BORRAR) 

DATATETCAZZA ICA ZABS TABS ZARZA 
LO ERAOO ERA AO DONNA 

PRINT 100,BORRAR 

FORMAT(ALO) 

END 

HEHE NN NENE INTI NINO NIN 


SUBROUTINE GOTOXY(I,J,PONER) 
CHARACTER*1 FILA(9) 
CHARACTER+*2 FILA1,FILAS 
CHARACTER*9 PONER 


100 


2010000060000 00n 


10 
100 


DATA FILA(1),FILA(2),FILA(3)/2*5C0”,2*1B”,2>5B”/ 
k=1=(1/10)*10—+—ICHAR(>0”>) 

FILA1=CHAR(I/10 + ICHAR(*0”*))//CHAR(K) 

K=J- (J/10)*10 + ICHAR(*0”) 

FILAS=CHAR(JI/10 + ICHAR(*0”))//CHAR(K) 
PONER=SFILA(1)//FILA(2)//FILA(3)//FILA1//”3*//FILA2//*H” 


RETURN 
END 


SUBROUTINE GOTO 

CHARACTER*1 FILACI1) 

CHARACTER+11 PONER 

EQUIVALENCE (FILA,PONER) 

A A y Y 

* 2*48*,2*00* ,2*00*/ 

PRINT-100,PONER 

FORMAT(1ih+,49) 

END 

SUBROUTINE —FRAME CUPPERLEFTX UPPERLEFTY,LOWERRIGHTI: LOWERRIGHY) 


PARAMETROS PARA MANEJAR EL” CURSOR: 
UPPERLEFTX + 
EJE IZQUIERDO X 


UPPERLEFTY: 
EJE TZ2QUIERDO Y 


LOWERRIGHX: 
EJE DERECHO PARTE DE ABAJO EN X 


LOWERRIGHY: 
EJE DERECHO PARTE DE ABAJO EN Y 


PONIENTE IET NNNININO 
INTEGER UPPERLEFTX UPPERLEFTY 

CHARACTER*1  LEFTCORNER,RIGHTCORNEL , LINEA, MURALLA 
CHARACTER*9 PONER 


CALL GOTOXY(UPPERLEFTX,UPPERLEFTY,PONER) 
LINEA=1th= 

LEFTCORNEL=1H 

RIGHTCORNEL=1H, 

MURALLA=1H 

PRINT 100,LEFTCORNEL 

DO-—10 T=UPPERCLEFTX+1,LOWERRIGHTX-1,1 
PRINT 100,LINEA 

CONTIMUE 

FORMAT (1H+,80A1) 

PRINT 100,RIGHTCORNEL 

DO-20 T=UPPERLEFTY+1 LOWERRIGHTY-1,1 
CALL GOTOXY(UPPERLEFTX,I,PONER) 
PRINT-100,MURALLA 

CALL GOTOXY(LOWERRIGHTX,I,PONE-) 
PRINT-100,MURALLA 

CONTINUE 

CALL GOTOXY(UFPPERLEFTX,LOWERRIGHTY ,PONER) 
LEFTCORNEC=1Ht 

RIGHTCORNEL=1HJ3 
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PRINT-100,LEFTCORNEC 
DO-—30 —T=UPPERLEFT+1LOWERRIGHTX—1,1 
PRINT 100,LINEA 
30 CONTINUE 
PRINT 100,RIGHTCORMNEL 
RETURN 
END 


SUBROUTINE LEERLINE 


ce 
(9) Esta rutina lee una linea de comando desde 
e el terminal TTY como máximo lee 80 
£ caracteres en LINEA. 
e 
COMMON —7COMMAND/LCIN 
CHARACTER+*80 LIN 
[e 
£ 


READ(%*, 100)LIN 
100 FORMAT(ASO) 
END 


SUBROUTINE SCANNER 


£ 
E Esta subrcoutina hace el scanner de la linea 
c introducida por la TTY luego se procede a 
€ efectuar el parse de la misma linea, 
C 
COMMON—7COMMAND/LC INEA, COMANDO 
COMMON /COMMAND1/LONGITUD, INDICE 
CHARACTER*13-—COMANDO 
CHARACTER+*1  TOKEN(13),LINEA(80),CR,TAB, BLANCO 
INTEGER CARACTER 
EQUIVALENCE-—(TOKEN, COMANDO) 
LOGICAL FLAG 
DATA CR,BLANCO,TAB/2”0D”,2*20*,2*09*/ 
a 
FLAG=.FALSE. 
LONGITUD=0 


DO 10 I=1,13 
10 TOKEN(T)=BLANCO 
e ignoro los blancos al comienzo, blancos y tabs 
I=INDICE 
J=1 
100 CONTINUE 
CARACTER=ICHAR(LINEA(I)) 
IF ((CARACTER.EQ.32) .OR. (CARACTER.EQ.9)) THEN 
I=1+1 
GO TO 100 
ENDIF 
200 CONTINUE 
IF (CICHAR(LINEACI)).6E.97).AND. CICHAR(LINEA(I)).LE.122)) 


THEN 
€ Convierto de minusculas a mayusculas y 
E ponga dende debe estar. 
TOKEN(J)=CHAR( ICHAR( LINEA(I) ) - 32) 
I=1+1 
i=j+t 
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000000000 00000000N00o 


GO TO 200 
ENDIF 
Verifico si es mayúscula. 


IF (CICHAR(LINEA(1)).GE.65).AND. (ICHAR(LINEA(I)).LE.90)) THEN 


Era mayúscula la guardo en lugar segurc. 
TOKEN(J)=LINEA(I) 
I=1+1 
J=J+1 
Go TO 200 
ENDIF 
No era letra mayúscula ni minuscula 
Veremos si es un número. 


FAA TECHARAL INEA CI 656,48) AND ¿LE HAROLINEA CIO LE-57))—PHEN 


Era un número la guardo en lugar segura. 
TOKENMCI)=LINEACI) 
I=I+1 
J=J+1 
60 TO 200 
ENDIF 

Veremcs que cosa extraña puede ser: 


IF ((CICHAR(LINEA(1)).E0.61).0R.(ICHAR(LINEA(I)).E0.42).0R. 


* CICHARICINEA CIO) E0.43).0R. CICHARILINEACI)).E0.47).0R. 
* (ICHAR(LINEA(I)).E£0.59) 
£.OR-CICHAR(CINEACIO).E0.43)) .AND-(J.EQ.1)) THEN 
Es-—un— signa *./ 1. too bien =: 
TOKEN(J)=LINEA(I) 


I=1+1 
INDICE=1 
CONGITUD=J 
REFURN 
ENDIF 

Era un delimitador. 

INDICE=1 

CONGITUD=I3=1 

RETURN 

END 


SUBROUTINE FILTROS(FICTRO) 


ESTA SUBROUTINA CALCULA DIFERENTES TIPOS DE 
FILTROS PARA EL MANEJO DE IMAGEN 

EL VALOR DE LA VARIABLE FILTRO NOS INDICA 
CUAL SERA EL FILTRO A GENERAR 


FILTRO=1 
ENTONCES SE CALCULA FILTRO BUTTERWORTH:. 


FICTRO=2 
ENTONCES SE CALCULA EL FILTRO EXPONENCIAL. 


FILTRO=3 
ENTONCES SE CALCULA EL FILTRO TRAPEZOIDAL. 


FICTRO=4 
ENTONCES SE CALCULA EL FIELTRO LAPLACIANO. 


AAA 


COMMON /BUFFER/ A 
COMFLEX-A(4096),AB(64,64) 


87 


200 


300 


400 


10 
500 


600 


FOO 


800 


35 


388 


EQUIVALENCE (A, AB) 
INTEGER+*2 FILTRO 
AQUI DEBE IR UN IF THEN ELSE. 


CALL CLRSER 

CALL GOTOXY(2,10,POMER) 

WRITE (+, 100) 

FORMAT(1H+,*"Deme el valor de las puntos de carte en ejes X e Y 


READ(*,200)CO0KTEX,CORTEY 

FORMAT(2F 10.0) 

IF (CORTEX.LE+0,0) CORTEX=1.0 

IF (CORTEY.LE.O.O0) CORTEY=1.0 

CALE-GOTOXY(2,10,PONER) 

WRITE (x*,300) 

FORMAT(1H+,*Deme el ancha en eje X de la transición 1”) 
READ(+*, 400) ANCHOX 

FORMAT(F10.0) 


FC ANCHOX LEO 0) ANCHOX=1,0 


CALL GOTOXY(2,10,PONER) 
WRITE (*,500) 
FORMAT(iH+, "Quiere filtro pasa alta (1) o filtro pasa bajos (e 


READ ,600)TIPO 

FORMAT(F10.0) 

TF A TIPOGNE 1070 CTIPO<NE-2.0)7)80-TO0-10 

CALL GOTOXY(2,10,PONER) 

WRITE(*,700) 

FORMAT (1H+,*Orden del filtro 2”) 

READ(*,800J)0RDEN 

FORMAT(F10.0) 

IF (ORDEN.LE.0.0) ORDEN=1.0 

DO-20—1=1533 

DO 20 J=1,33 

11=32+1 

12=34-1 

J1=32+J3 

J2=34-3 

DX=FLOAT ( I-1)/CORTEX 

DY=FLDAT(I=1)/CORTEY 

D=CABStEMPLXCDX,DY)> 

IF (FILTRO.NE.3) THEN 
IF(ORDEN.NE.1.0.AND.D.NE.0.0) D=D**ORDEN 
x=0.0 
IF (FILTRO.EO0.1) X=1.0/(1.0 + 0.412xD) 
TF-0FILTRE0.2— AND. D ET 100) X=EXP(=0.347 
*D) 
TF FICTRO:E0:4-)>X=D 
ENDIF 

X=ORDEN 

1F<D7EE-150)—60-—0-2e 

D=(D-1 .O)*CORTEX 

X=1.0——— ORDEN 

IF(D.GT.ANCHOX) GO TO 22 

X= (ANCHOX=D)7ANCHOX 

IF(TIPO.E0.1.0) X=1.0 — X 

IF (ORDEN.NE.1.0)X=X*EXORDEN 

CONTINUE 

IF (D.EQ.0.O0 .AND. FILTRO.NE.4) X=TIPO 


£ LLENO” AHORA COS PIXEL”S 


IF (J1.E0.65 ) GO TO 40 
TEC IES 6700 TO 41 
AB(T1,J1)=X 


41 IF (12.E0.65) GO TO 40 
AB(12,3J1)=X 
40 IF(-J2.E0.65-)-G0-TO-20 


IF( 11.E0.65 ) GO TO 42 
ABI, J2)=X 


42 IF(12.£0.65) GO TO 20 
AB(12,32)=X 

20 CONTINUE 
RETURN 
END 


SUBROUTINE TRANSFORCINDICE) 
COMMON 7/BUFFER/—A 
COMPLEX A(564,654),A41(40965) ¿WORK (64) 
EQUIVALENCE CA,AR) 
REAL AAA(564,64) ,RA(4095) BE (64,64) 
EQUIVALENCE (A,AA,AAA,RA), (AAC2049),BB) 
1F-CINDICE.EO.1)—THEN 
e TRANSFORMADA DIRECTA DE FFT. 
PRINT*,* Transformada directa de fourier (FFT)3” 
SIGNO=+1.0 
K=1 
DO-—100—T=1t,64 
K=3-K 
DO-100—J=kK,64,2 
100 A(1,J)=-A(1,J) 
DO-110 1=1,64 
CALL FASTFFT(A(1,1), SIGNO) 
110 CONTINUE 
DO—120-—J=1,64 
DO 125 I=1,64 
WORKCI)=ACI, ID) 
125 CONTINUE 
CALL FASTFFTIWORK, SIGNO) 
DO 126 I=1,64 
ACI, I)=WORK(I) 


126 CONTINUE 

120 CONTINUE 
k=1 
DO 130 I=1,64 
K=3-K 


DO 130 J=k,64,8 
ACIJIJ==A(1,3) 
130 CONTINUE 
RETURN 
ELSEIF (INDICE.E0.2) THEN 
= TRANSFORMADA INVERSA DE FOURIER. 
SIGNO==1.0 
k=1 
. DO 200 I=1,64 
K=3=K 
DO 200 J=2,64,2 
200 AT,I)==A(1,3) 
DO 210 1=1,64 


223 


226 
220 


230 


101 
102 


103 
100 


130 
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CALECFASTEFTCACI,TIT),SIGNO) 
CONTINUE 

DO-220-J=1,64 

DO 225 1=1,64 
WORKE(I)=A(J3,1) 

CONTIMUE 

CALL FASTFFT(WORK, SIGNO) 
DO 225 I=1.54 

ALT, 1)=WORK CI) 

CONTINUE 

CONTINUE 

K=1 

DO-230-I=1,64 

K=3-k 

DO-230-J=K,64,2 
ACI.J)=-A(1,J) 

CONTINUE 

RETURN 

ENDIF 

END 


SUBROUTINE FASTFFT(A, SIGNO) 


COMPLEX—A(64),U,DOBLEW,TE 
DATA NM1,NV2,N,M,P1/63,32,64,6,3.141592653589793/ 
J=1 

DO 100 I=1,NM1 

TF-CGE-3)-G0-T0-101 

TE=ACI) 

ALI)=A(I) 

ACII=TE 

E=NV2 


IF US GE- 3) 60-T0-103 
JI=IK 

E=EZ2 

GO TO 102 

I=IFK 

CONTINUE 

DO 110 L=1,M 
LE=S2x%x*L 

LE1i=LE/2 

U=(1.0,0.0) 
ANGULO=PI/FLOAT(LE1) 
DOBLEW=CMPLX (COS (ANGULO) ,-SIGNO*SIN(ANGULO) ) 
DO 120 J=1LE1 

DO 115 1=J,N,LE 
IP=I+LEI 

TE=A(1P)*U 
ACIP)=A(I)-TE 
ACII=ACIATE 
CONTINUE 

U=UXDOBLEW 

CONTINUE 

CONTINUE 
DO-130-1=1,64 
ACI)=A(1)/8.0 
CONTIMUE 

END 
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SUBROUTINE—INDUTCOPCIONN BUFFER, CANAC) 
COMPLEX BUFFER(N) 
INTEGER OFCION,CANAC 


LARGO=54%N 
IF COPETIONEO 1) THEN 


READICANAD) CBUFFERCEO)K=kE17K2) 
CONTINUE 

RETURN 

ELSEIF (OPCION.EO.2) THEN 

DO 200 I=1,LARGO, 64 

K1=1 

k2=1+63 

URITE(CANAL) (BUFFERIK), K=KI¿K2) 
CONTINUE 

RETURN 

ENDIF 


FPRINT *," ERROR: Na existe la cpcian que me”, 
x indico en uno de lus parametras” 
END 
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NIVELES MIN, MAX 3 —0:;0000000E+00 1,.000000 


-32 
-31 
=30 
-29 
=28 
-2? 
=26 
-25 
=24 
-23 
=22 
-21 
=20 
-19 
-18 
=17 
-16 
=15 
-14 
== 
=12 
=H1 
-10 

-9 

-8 

=7 

-6 

-=5 

-4 

=3 

=2 

= 


Pbhebn 
VN-Po00b0JYJNIDADmurOo 
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TAILS AAA AAA 
ARALL LL LLL LLL LLL ALL L LLL LLL ALLE LLL LLL LLL LLL LLL LLL LLL LL LL LLL LLL 
ALA AAA LA ELA A LALALA LALALA 
BAALA LL LLL ALL LLL LLL ELL LLL ALL ALL LLL LL LLL ALL LLL LLL LL LL LLL LLL LLiLiL 
IA AAA AA AAA 


AAA A AAA AAA AAA LLL 
AAA LL LLL LLL LLL LLL LLL LLL LLL LLL LL ELL LLL LLL LLL LL LL LLL LL LL LLiLit££ 
ALLA LALALA AAA LALALA ALA ALLA A LL LLL LAA 
ALALL LLL LLL LL LLL LLL LLL LLL LALA LLL LLL LLL LL LL LLL LLL LLL LLL LL LLL 
ALLA AAA LLL AAA LLL ALL LLL ALLA ALLA LLL 


ELLLLL LLL LLL LLL LLL LLL LLL LLL EL LLL LLL LLL LLL LLL LLL LLL LLL LLL LLL LLL 
AAA LA LL LLL LLL LL LLL LLL ALL LL ALL LLL LLL LLL LLL LLL Li LL LALA 
ALELLA LLL LLL LLL LLL LLL LLL LLL LLL LLL LLL LIL LLL LLL LIL LIL LLL LLLILLL 
AAA AAA LALALA LLL LL 
ALLE LLL LLL LLL LLL LLL LLL LLL LLL LLL LLL LLL LLL LLL LLL LILA LLL LLL LL LLL 


ALLA LL LLL LLL LLL LLL LLL EL LLL La LLL LLL LLL LLL LL LLL LLLiri si Litinit£ 
LALA AAA AAA A ALLA A AA 
ALELLA LLL LLL LE LLL ELLE LLL LLL LLL LL LLL LLL LLL A LLL LLLLiLto 
AAA A ALL LA LALA LLL 
AAA ALL LLL LL LL ELL LL LLL LL LLL LL LLL LLL ALL LL LLL AL LL Lab iLi££it 


ALLA LLL LALA LLL LLL LLL LLL LL LL LLL AL LA LL AL LALALA A ALLA LLL LLL LA 
ALLA LALA LLL LLL A ALLA LAA LALA LLL LAA LLL 
AAA LA LALALA LLL LLL LLL LLL LALALA LLL ALL LA LLL LLL LLL ALA LLLLL 
ALLA LLL LLL LLL LLL LLL LLL LLL LLL LLL LLL LLL LLL LLL LLL LLL LALA LL LLL 
AAA AAA AAA AAA AAA AAA AAA AAA AA LLL La 


LLL LALALA ALLA LALALA LLL LLL 
ALLA ALL LL LLL LLL LLL LL ALL LL EL LLL LL LLL aL LL LLL LLL LLL ari iLiLiiiii£t£ 
ERA ALA AAA LALA LLL 
ALLLLLL LLL LL LLL LLL LLL LLL LL AL LLL LLL LLL LLL LLL LLL LLL LL LLL iLiLititit 
AAA LALALA LL ALL LLL ALL LALA LLL LLL 


ERA LLL ALLA LL LA LLL LALA LLL ALL LLL AL 
ALLA LLL L LLL LLL LLL LLL LL ALL LLL LLL LLL ALL LLL LLL LLL LL LLL LL otitit 
ALLE LLL LLL LLL LLL LAA LLL LLL LL LLL LLL ALLI ILLA 
TASAS SAA AAA AAA AAA AAA ALL 


Fig. 1. Test de patrones horizontales. 


NIVELES MIN, MAX 3  0,D000000E+00  1156.000 
=32 
-21 
=30 
-29 
=28 
-27 


Fig. 2. Módulo de la transformada de Fourier de la Fig. 1. 
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NIVELES-MINS—MAX 70 D000000E+00 


—3e 
== 


=35 


Fig. 3. Test rectangular. 


1.000000 


LLL 
LLLLLLLLLL 
LALALA ALLA 
ALLLLLILLL 
LLL LALALA 
ILLLLLLLLL 
LLL 
LLLALLLLLL 
AA AAA 
LLLLLLLLLL 


NIVELES MIN, MAX: —0.0000000E+00—2.441406 
=32 
-21 
=20 
-29 
=28 
-27 
=26 
=25 
-24 
-23 
-22 
=21 
-20 
=13 
-13 
=17 
-16 
=15 


-4 + 
== HUA 
-2 FUESE ESUA 
=> FHHESLIEASE 
HESLLL EIA 
HEELIEAT 
HUB ESIA 
HIIARAA 
++..+ 


31 
Fig. 4. Módulo de la transformada de Fourier de la Fig. 3. 
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APENDICE 4 


GLOSARIO 


Imagen. Representación espacial de un objeto, escena u otra imagen. 
El nivel de gris I(x,y) correspondiente a la coordenada espacial (x,y) es, en 
general, proporcional a la energía radiante. 

Nivel de gris. Corresponde al número o valor asignado a la posición 
(x,y) de una imagen. Este número es proporcional a la reflectancia o trans- 
mitancia de una pequeña área denominada pixel. 

Intensidad radiante. De un punto o un objeto es una medida de la 
potencia radiante por estereorradián radiado o reflejado por un objeto. La 
intensidad radiante es una función de la naturaleza del objeto, ángulo de 
observación, longitud de onda espectral y ancho de banda. 

Reflectancia o coeficiente de reflexión. Es el cociente entre la ener- 
gía por unidad de tiempo y por unidad de área (densidad de potencia ra- 
diante) reflejada por un objeto y la energía por unidad de tiempo y por uni- 
dad de área incidente en el objeto. 

Transmitancia o coeficiente de transmitancia. Es el cociente entre 
la energía por unidad de tiempo y unidad de área (densidad de potencia 
radiante) transmitida por un objeto y la energía por unidad de tiempo y uni- 
dad de área incidente en el objeto. 

Contraste. El contraste de un punto frente a un fondo puede medirse 
por: 1) su razón de contraste, que es el cociente entre el valor mayor de 
la transmitancia del objeto o del fondo y el valor menor de la transmitan- 
cia del objeto o del fondo, y 2) diferencia de contraste, que es la diferencia 
entre el valor mayor de la transmitancia del objeto o del fondo y el valor 
menor de la transmitancia del objeto o del fondo. 

Densidad. La densidad en un punto (x,y) de una fotografía es una me- 
dida de la capacidad de absorción de la luz de la plata o del tinte deposi- 
tado en cada posición. 

Densitómetro. Dispositivo utilizado para medir el valor medio de la 
densidad de la imagen correspondiente a una pequeña área de una foto- 
grafía o de una transparencia. Cuando el área es menor que unos cente- 
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nares de micras cuadradas, el dispositivo recibe el nombre de micro- 
densitómetro. 

Resolución. El término resolución hace referencia a cómo un siste- 
ma, proceso, componente o material puede reproducir un objeto aislado 
o discriminar entre objetos próximos. El límite de resolución o resolución 
espacial viene dada por la dimensión más pequeña del objeto que puede 
ser discriminada. 

Poder de resolución. El poder de resolución de un sistema de trata- 
miento de imágenes es una medida de su capacidad de discriminación en- 
tre objetos próximos. El poder de resolución se mide en pares de líneas 
por milímetro, es decir, el máximo número de líneas y espacios por milí- 
metro que pueden reconocerse. 

Imagen de vídeo. Es una imagen que presenta un formato capaz de 
ser presentada en un tubo de rayos catódicos. 

Cuantificación. Proceso mediante el cual cada nivel de gris de una 
imagen, fotografía, etc., se le asigna un valor entre un conjunto finito de 
niveles de gris. 

Filtro espacial. Transformación de una imagen utilizada para redu- 
cir el ruido o mejorar algunas características de la imagen. 
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El procesado de imágenes es un campo de 
reciente y rápido desarrollo, con 
importantes aplicaciones en áreas tan 
diversas como la mejora de imágenes 
biomédicas, robóticas, teledetección y 
otras aplicaciones industriales y militares. 
Se presentan los principios básicos, los 
sistemas y las técnicas de procesado más 
usuales. 


